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THE INFLUENCE OF THE INTERNAL STRESS IN 
DIELECTRICS ON THE POLARIZATION OF 
REFLECTED LIGHT 


By Zygmunt (Ce esol) OA Institute of Theoretical P 
D of Wrocław 


(received July 7, 1950) 


The refractive index considered as a function of the strains has been introduced 
into the equations decribing the state of polarization of light reflected from a strained 
dielectric. Glass has been chosen as example since the Neumann constants are 
known from experiment for this material only. The relation between the rotation 
of the plane of polarization of reflected (plane polarized) light and uniform traction 
normal to the plane of incidence has been deduced, as-well as a relation between 
the rotation of the plane of polarization of reflected light and the state of plane 
stress within the dielectric. The latter relation yields one equation for the determi- 
nation of the two principal stresses. The second equation needed can be deduced 
from the fundamental equation of photoelasticity connecting the phase shift of trans- 
mitted light with the stress within the dielectric. 


For transmitted light the influence of pressure on polarization 
is known for a long time. About 1815 Brewster observed that optically 
isotropic bodies under one-sided pressure behave like negative uniaxial 
crystals, and under one-sided tension like positive uniaxial crystals. 
He found also the phase difference between the ordinary and the 
extraordinary rays to be proportional to the traction. In 1841 F. Neu- 
mann published a theoretical investigation of the above problem based 
on the elastic theory of light and leading to the following formulae 
connecting the velocities of light in the principal directions of a strained 
medium with the principal displacements: 


V1 = Vo t+ Plot e3) +941 

Vg= Vt ples t 1) +4 ee _{1) 

Pas vot P(t e2) +463. i 
Here v; is the velocity of a waye polarized in a plane perpendicular 
to the 7-axis, p, q are s N eumann’s constants, and e, €z, €g are the prin- 
cipal elongations. 
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The influence of traction on the polarization of reflected light 
was found experimentally in the region of complete reflection by Volke 
(1910), and in the vicinity of Brewster’s angle by Lummer and Sorge 
(1910). 

Recently problems connected with the subject of the present 
paper excite again more intense interest due to the development of 
technical photoelasticity, based on optical methods of measurements 
of interior stresses. This may be illustrated, for instance, by the fact 
that the list of papers and books on elastooptics published between 
1935—1946, contained in Frocht’s text-book on Photoelasticity (1948) 
and supplementing the bibliography given in the fundamental Treatise 
on Photoelasticity of Coker and Filon (1931, p. 190), comprises 14 full 
pages. 

It seems, however, that the technical methods used now in photo- 
elasticity require some essential supplements. The present paper is 
an attempt to introduce such a supplementary method based on the 
investigation of reflected light. 


Relation between the polarization of reflected light and 
uniform traction exerted on a plate 


Consider a plate of dielectric in the form of a rectangular parallel- 
epiped. Under uniform pressure parallel to the z-axis, the plate will 
become an anisotropic body having optical properties of a negative 
uniaxial crystal. The direction of the pressure P corresponds to the 
optical axis. We apply the general equations describing light pro- 
_pagation in anisotropic media (Schulz, 1926, p. 685) and assume the 
‘incident light to be perpendicular to the direction of the pressure. 
The principal axes will be designated by the indices 1, 2,3. Taking 
into account the axial symmetry of our problem, the fact that our 
system is referred to principal axes, and that the uniradial azimuths 
are 0,=0, 0,=2/2, we obtain the following relations between the 
amplitudes of the incident and the reflected light 

RaR 22am) . (peter OF (2) 

n (p+ x3) 

where is the angle of incidence, y,, ys are the angles of refraction, 
E., E the amplitudes of the electric vector of the incident light 
Tee Rs the amplitudes of the electric vector of the refracted light. 
Index 1 refers to the components in the direction of axis 1, index 3 
to the components in the direction of axis 3. For the ordinary ray 
the electric vector oscillates in the direction of axis 1, for the extra- 


=~ j 
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ordinary one in the direction of axis,3. Since the glass plate assumes 
the properties of a uniaxial crystal „v <v and n,>n3, where vı and v 
are the velocities and n, and n the refractive indices of the ordinary 
and the extraordinary rays respectively. Owing to -axial Symmetry 
around a direction perpendicular to the optical axis, i.e., to the di- 
rection of the traction, the intersection of the plane of incidence with 
the wave surface consists of two circles. The direction of wave pro- 
pagation is the same as the direction of energy propagation. The in- 
cident and refracted rays are in the same plane as the normal to the 
refracting surface and the index of refraction of the extraordinary 
ray does not depend on the angle of incidence. Thus, in our case, the 
two usual laws of refraction are also applicable to the extraordinary ray. 
Since . 


= <<1, where Ajw= Li Vo, 


we can write 


Ae IAW aie c rt n =n — 2) 
TU o Ayw ARN Aw vr 
Vo 


and hence from Neumann’s equations (1), we have 


nı =n[1—p (e, + e3) — g'e], (3) 
ny=n[1— p (e+ ey) — q'es], 


where 


In our case 
EE TO (4) 
4] a 
when 7 is the elasticity modulus and u Poisson’s coefficient, and equa- 


tions (3) can be written in the form 


ni=n+O,P (i=1,3), rion 
where 

= 7 (= Hp" + MAIN ar 
i. (6) 
C3= ple Sey pil 


It has been already stated that both the ordinary and extra- 
ordinary rays comply with the usual laws of refraction, we have 
i . a 
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therefore 


guy _, — [msing G13): (7) 


sin Kin. COS Xi 1; 


Since in (5) O,P<<1 we can write 


mi=n?+2n0,P (t=1,3) (8) 


Vn sinp = =N 14 A k (9) 


and hence 


where we have put 
n?— sip =N. 


Now we can express the ratio R,/R, as a function of the pres- 
sure P and the angle of incidence p, as follows. From (2) we have 


kı E, sin (p—ya) 608 (p+ %1) sin (P+ Xs) 


PEO (G77) sin (p + 71) sin (p— val 


Inserting here the values of sin y, and cos y, from (7) and taking 
into account (5) and (9), we have 


(senate eer P42) a 
eae e 


Finally, neglecting products of the small quantities nC,P/N and 
nO0,P/N, we get 


Rene er 
R= 0, GP EH’ (10) 


where 


C—O. 
(cos g Vn? —sin? pg —sin? p |a (0+0) + fate cos o| 


C 
a 
Vn?—sin? y 
F= (n?—1) (cos g Vn?—sin? g —sin? p), 


G= (COS o n?— sin? pọ + sin? p| 0+0) n COs o| 
—sin? y | 


C 
Fna i a co 
; Vn? —sin? y 2 
H= (n?—1) (cos pVn?—sin? y + sin? p) 
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By way of example the two following tables!) have been dressed 
for normal incidence, y=0. In these tables 
E R 
a tana = =o 
E,’ 0 n R, (11) 
Table I. 
For np =1,50750, n= 4800 kg/mm?, u=0,274, p’ =0,274, q’= 0,166 


C,=0,0000482, C,=0,00000497 ` 
Ee SR | CP 


175° 0’ 0” 1,50750 1,50750 
174° 59’ 58’| 1,50760 1,50751 
135° 0’ 0” 1,50750 1,50750 
134° 59’ 48”| 1,50760 1,50751 


Table II. 


For np =1,96250, ņn=5035 kg/mm?, u=0,26, p’=0,427, q’=0466 


C,=0,0000759, C,=0,0000951 


Ny 


175° 0’ 0” 1,96250 1,96250 
175° 0 1” 1,96265 1,96269 
135° 0’ 0” 1,96250 1,96250 
135° 0’ 3” 1,96265 1,96269 


If we compare this with the anisotropy of such crystals as 


Iceland spar: n,=1,6530, nz=1,4841, 

quartz: N,=1,5410, n, =1,5500, 
we must conclude that the anisotropy of glass is extremely week even 
under very high pressures. It may be surmised that for some other 
photoelastic materials (see below) this anisotropy might be much greater 
than for glass, though the Neumann constants of these materials are 
as yet undetermined, as far as it is known to the author. 


Polarization of reflected light in the case of a plane stress 


In photoelasticity the state of plane stress is treated more or less 
as follows. At any point there exist two mutually perpendicular di- 
rections, called principal directions of stress, in which the stress is 


1) Values of constants have been taken from a paper of Pockels (1902) and 
refer to two sorts of optical glass: S 205 and S 57 of Schott respectively. 
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a maximum or a minimum. The transmitted- light is investigated in 
a rectangular coordinate system with axes pointing in these directions. 
A beam of plane polarized light is normally incident on the surface 
of the plate. After passing through the plate it becomes elliptically 
polarized. The phase difference and the orientation of the axes of the 
ellipse (and thus also of the principal directions of stress) can be found 
experimentally. 

The formulae for the phase difference may be found from Neu- 
mann’s equations (3). In the case of a state of plane stress, one of the 
principal stresses, say Pa, vanishes and 


i 1 
=P. + HPs) a=" (Pit Ps), Ce Es ad (12) 


(see, e. g., Love 1906, p. 204, where the axes of coordinates have been 
taken in opposite directions, which explains the difference in sign of 
the right-hand membres. The same applies to our formulae (4), as 
a special case of (12)). Inserting these values in (3) we get 


with the same values of the contains ©, and C; as in (5). Hence 
Ny —Ng = (C3 —C,)(P,—P3) (14) 


and the phase difference is given by 
d 
A= 7 (03—01) (Py —Ps) = k(P, —Ps). 2 (15) 


This is the fundamental equation of photoelasticity; d is the thickness 
of the plate, 2 the wave-length, and 


3 (0—0): (16) 


If we express k in mm?/kg 
k*¥1 for bakelite, phenolite, 


0,2 „ celluloid, 
0,05 ,, glass. 


Thus we see that the value of k for glass is very disadvantageous; 
at equal stresses the phase difference, and consequently also the dif- ; 
ference of the refractive indices, is for glass two order of magnitude 
smaller than, e. g., for phenolite. This means that in the case for which 
the calculation has been carried out in the preceding section the- 
rotation of the plane of polarization should be much more psroaptibla 
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for phenolite. Unfortunately, as already mentioned, an exact com- 
putation cannot be carried out since Neumann’s constants have been 
determined, as it seems, only for different sorts of glass. 

Equation (15) is Ane one equation for two unknowns. The second 
equation can be obtained by different, rather complicated, optical or 
mechanical methods, which are said, however, to be comparatively 
inaccurate (Coker ane Filon 1931, Maa 1939, Frocht 1948). 

According to the work in the preceding igor the second equa- 
tion can be obtained also by investigating the A aon of the 
refleçted light. In our case, for normal incidence and for a material 
in a state of plane Stress, Fresnel’s formulae read 

Ny— AEN al 
Lh nyo ng +1 

The amplitude of the incident light is now resolved not with 
respect to the plane of incidence (which in the case of normal inci- 
dence is only a limiting notion) but with respect to the principal axes 
of stress, which are in general different at different points and can 
be found experimentally by investigating the transmitted light .by 
means of a compensator and an analyser. 

Inserting the expressions (13) for the refractive indices and ne- 
glecting the small products C,P,, CPs, C;P, and 0}P}, we get from (17) 

Ry, E (aP, +bP}+ m) 


ie ee he BL 18 
i. E (bP + aP} +m)’ aS) 


On ee 


F}. (17) 


where 
a=(n—1)C,+(n+1)C3, b=(n+1)0,+(n—1)03,, m=n?—1. 


Making use of the notation introduced in (11) we can therefore 
write 
AP,+BP,+C=0, (19) 
where 
A=atana+btane, B=atane+btana, C=tana+tano. (20) 


A,B,C are given functions of known, or immediately measurable 

‘quantities, 

| Thus equation (19) together with equation (15) form a set of 
|two linear equations from which the two principal stresses P, and Ps 

‘may be calculated, unless the determinant —k (A+B) of these equa- 
tions vanishes. This is the case either when k=0 or A+B=0. From 
(20) we see that A——B implies either tan a=— tan ọ or a=—b. 

[In the former case C=0, P,=P, and, by (15), the phase difference A 
‘vanishes. The latter case PRA, Ops suc and from (13) we get then 
,+n,=2n. This means that there would be no rotation of the plane 


a 


J 


( 
ie 
! 
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of polarization due to stress in a body whose refractive index under 
no stress would be equal to the arithmetic mean of both the refractive 
indices under stress, irrespective of the magnitude of the stresses. For 
Neumann’s constants we would obtain in this case from (6) the relation 


(3u —1)p =4(1— u). 


It seems, however, that the case C= — 0; is never realized in nature 
since hitherto both coefficients C, and 0, have been always found to 
be positive. 

The determinant vanishes also for k=0, i. e., for C,=C,. Accor- 
ding to Pockels (1902) this occurs, e. g., for glass with a refractive 
index of n=1,88, but this is a quite exceptional case. 

Thus we see that our equations always have a solution except 
when we do not find any phase difference of the transmitted light or 
any rotation of the plane of polarization of the reflected light. 

Finally, it may be noticed that altough at present the method 
proposed here can give only inaccurate results due to the impossibility 
of attaining a sufficiently high degree of anisotropy, yet the rapid 
progress in the production of plastics for use in photoelasticity entitles 
us to hope that this method may prove useful in the future. 

My thanks are due to Dr R. Ingarden for suggesting the problem 
and for many helpful discussions. 
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THE LOWERING OF THE FREEZING POINT IN THE 
THEORY OF DIPOLAR COUPLING 


By A. PIEKARA, Physical Laboratory of the Institute of Technology, 
Gdansk 


(received July 10, 1950) 


The theory of dipolar coupling developed by the author in a previous paper 
(Acta Phys. Polonica, 10, 107 (1950)) has been applied to the calculation of the 
lowering of the freezing point of dilute solutions of polar substances in non-polar 
solvents. The lowering of the freezing point in these solutions (AT) is less than 
that in „ideal“ solutions (A7ja) in which the dissolved molecules do not exert any 
forces on each other. The quantity j=(4 T—ATia)/ATia has been calculated on the 
assumption of a ,,coupling of the second kind“ tending to associate the dissolved 
molecules in antiparallel couples (for details, see the above mentioned paper). The 
value of (dj/dm),,-9, where m is the molar concentration, determines the tangent 
to the curve j(m) at m=0; it is found to depend on the shortest possible distance ro 
between the molecules. 

The theoretically calculated tangents fit well enough with the experimental 
curves. It must be emphasized that for nitrobenzene, orthonitrophenol and para- 
nitrophenol no new constants have been introduced, as the value of r,=4,6 A has 
been taken from the above mentioned investigation, where it was determined from 
a discussion of the dielectric, electro-optical and magneto-optical phenomena in 
solutions of nitrobenzene in non-polar solvents. 


1. Introduction 


It was found in a previous investigation (Piekara 1950) that 
three kinds of dipole couplings are practically sufficient to account 
for the behaviour of the four phenomena of molecular orientation in 
solutions of polar substances in non-polar solvents. Notwithstanding 
_ the large differences in nature between these phenomena (for instance 
between the dielectric polarization and the magneto-optical _birefrin- 
gence) a single group Of constants is sufficient to describe them all 
in conformity with the experimental data. The purpose of this paper 
_is to apply the aforesaid theory to the explanation of a phenomenon, 
in which the molecular orientation does not actually play any part 
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whatever, but which is nevertheless subjected to the action of inter- 
molecular forces, namely the lowering of the freezing point of very 
dilute solutions of polar substances in non-polar solvents. 

Let us imagine that the molecules of the dissolved substance do 
not exert any force on each other and are, therefore, completely free. 
In such a case the solution is said to be ,,ideal”, and the lowering of 
its freezing point is the greatest possible. This lowering is expressed 
by the formula 

AT ia= îm, (1) 


where 2 can be calculated thermodynamically and, when benzene is 
the solvent, amounts to A=5,091. The molar concentration m is de- 
termined as the number of moles of the dissolved substance in 1000 g 
of the solvent. 

Let us imagine now that forces of cohesion develop between the 
molecules and the solution ceases to be ideal. Then, the osmotic pres- 
sure in the solution will decrease from Pig to P and involve a cor- 
responding decrease of the amount by which the freezing point is 
lowered (from AT to AT) is in volved. For an infinitely dilute 
solution the following relation following from thermodynamical consi- 
derations may be used: 

. AT ia AT B ia — P 
aT (Oa E RA 


j is a measure of the deviation of the given solution from the ideal 
state and, at the same time, it may be considered as a sort of meas- 
ure of the coupling forces between the molecules of the dissolved 
substance. The curve j(m) can be drawn fromex perimental data and, 
on the other hand, aalculated theoretically. In the present pa- 
per only very diluted solutions are to be considered and therefore 
only the part of the curve near the point m=0, i.e. the tangent 


(2) -© , will concern us. 
dm m=0 


iS 


2. The Calculation of Gal 5 
l dm m=0 : 
The problem of the calculation of j(m) is thus reduced to the 
calculation of the osmotic pressure P in the real solution, and is similar — 
to the calculation of the correction of the pressure in the equation 
of van der Waals (Herzfeld 1925, pp. 167 and 162; Fuoss 1934). i 


1 This is the mean of the values 5,065, 5,075 and 5,13 (see Tables Annulles 
des Constantes et Données Numeriques, 1941). ; 
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The free energy f of the molecule considered as a point particle 
in a very rarefied gas or a very dilute solution, is given by the formula 


BET pe 
eor) | Hav, (3) 


where ọ(T) is a function of temperature (with the exact form of which 
we are not concerned here), U is the potential energy of the molecule, 
and the integration is extended to all the positions of the given mo- 
lecule in the volume FV. 

It has been already assumed in the previous paper that the 
energy U consists of three parts, namely: 

(1) The energy of the coupling of the first kind; in dilute solu- 
tions this is the coupling of the dissolved molecules with the sol- 
vent. Since this part of the energy is independent of the coordinates 
of the molecule, the expression including it in the exponent can be 
put outside of the integral sign, thus causing only a change in the 
function g(Z), which does not influence the osmotic pressure. 

(2) The coupling energy of the second kind; this coupling acts 
between the given polar molecule and its nearest dipole neighbour. 
It has a strong influence on the osmotic pressure. 

(3) The coupling energy of the third kind, which consists of the 
interaction energy of more than two molecules and can be observed 
only for higher concentrations. Such cases will not be considered here 
since a very dilute solution is assumed. 

Therefore in formula (3) U indicates the coupling energy of the 
second kind only. This energy, being the energy of interaction of two 
dipoles, diminishes as the cube of the distance between the dipoles. 
Then it can be assumed that for a certain distance r>R the energy U 
does not differ materially from zero. Only more ,,intimate coupling” 
within the spheres of action of radius F will be therefore taken into 

account; by spheres of action are meant spheres enveloping the mole- 
cules, the coupling outside of them being disregarded. 

: The calculation of the free energy will be worked out only for 
‘infinitely dilute solutions. It can be assumed therefore that the mean 
distance between the molecular dipoles is very large compared with 
the radius of action R. Under such conditions only very few molecules 
are involved in more intimate coupling, and the grouping of molecules — 
by threes would be quite exceptional; this case will therefore be ex- 
eluded from the calculation. 

Let us assume that the solution is produced by adding the mo- 

ecules consecutively to the solvent. After adding the n-th molecule 
ne space being to its disposition will consist of two parts: the po- 
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4 A 
tentialless region of volume v—n) 5 R? outside the spheres of 


action and (n—1) spheres of action, inside of which TO 

The parts of the spheres which overlap are small and can be 
omitted. 

Consequently the integral of formula (3) for this molecule takes 


the form 
U 


V—(n—1) 2 B+ (n—1) |e PAV =V + (nV), (4) 
l r<R 
where J is the integral 
U 
J= f (iay (5) 


r<R 


taken over the sphere within which interaction is assumed to occur. 
The free'energy of the n-th molecule, omitting the expression depen- 
ding only on temperature, will be therefore: 


j= ETMV + =J= kT al 7(1+ =a). (6) 


n—1 
wale 
comparison to 1, it is therefore possible to restrict oneself to the first 
term of the development 


For very diluted solutions the expression J is very small in 


fo—kD TE J. (7) 


The expression for the free energy of the whole system will be ob- 
tained by adding up the foregoing expressions for all the molecules 
from n=1 to n=n together. Hence 


(n—1)n 


F=—nkT lnv— OV kTJ. (8) 
: OB ire: 
The osmotic pressure E= ae will be therefore 
nkT n 
where n has been substituted for »—1. Since TERAN the value 


V 
of j may be calculated from (2): 
nN 


j=57J- (10) 
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For an infinitely dilute solution the following relation is applied 


n  Nmd 
V 1000’ 


where N is Avogadro’s number and d is the density of the solvent. 
Thus we obtain the desired expression for the tangent to the curve j(m) 
at the point m= Q: 


dj Nad 
Cabs 0000" (1!) 


This result is identical with the one given by Fuoss (1934). The 
difference between the method of Fuoss and the present one will 
appear in the estimation of the integral J, for which it will be neces- 
sary to introduce some simplifying assumptions. The assumptions of 
Fuoss restricts the application of his theory either to solutions of 
electrolytes in non-polar solvents, or to a case when the molecules 
of the dissolved substance are spherical (Fuoss 1936). On the other 
hand the present author applies the same assumptions as in his pre- 
vious work (l. c.) concerning the phenomena of molecular orientation, 
namely, that the molecular dipoles are subjected to coupling of the 
second kind tending to an anti-parallel position. This theory can be 
applied, therefore, to all solutions of polar substances whose curves 
of dielectric polarization versus concentration are similar to that of 
nitrobenzene (the curve going down). With slight modifications it could. 
be applied also to solutions of dipole substances which show coupling 
of the second kind of a nearly-parallel type. 

In accordance with the foregoing, we put in formula (5) for nearly- 
antiparallel coupling 


2 
U=—5 cos 0, (12) 
and therefore 
J= fer 1)av, (13) 
. r<R 
where 
ss (14) 
| OS ap 
In a different form 
. R m P 
— cos O T, 
Fate f ryar fe —1) sin 6 dé, (15) - 
ro 0 See 
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where r, indicates the shortest distance between the coupled dipoles. 
The second integral can be calculated at once, giving 


R 
7, a 
J= x f (5 sinh —1) rar: (16) 
a i 
This integral may be calculated graphically, and then from formula (11) 


(2) follows. 
dm m=0 


3. Comparison with experimental results 


First we deal as in the preceding work with the solutions of 
nitrobenzene in benzene. We put: u= 4,23.10—!8 and T= 278,69 K (the 
freezing point of pure benzene and fie also of infinitely diluted 
solutions of other substances in benzene), or a= 0,466.10—%. The plot 
of the function under the sign of integration in (16) is shown in Fig. 1. 
The curve is going down so fast that the value of the integral J is 


Calculation of 


e. 
J- off sinh, -4r'dr, 


A= 0,466 10%” 


100 r= 46 A 


50 


4,95 “6l07 8 DG 10 11,42 773 r{A) x 


Fig. 1 


practically independent of the choice of the upper limit; it depends 
on the contrary to a very large extent on the value of the 1a limit 
so that a very slight change of the value of r) can produce apres vie 
or disagreement of the theory with experiment. The value of r = 4.6 A 
for nitrobenzene has been taken from the foregoing work; it haa pro 
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ved to be the best to account for the different phenomena of mole 
cular orientation, as the dielectric constant, the electric saturation 
and the electric and magnetic birefrigence. The value J = 2,50.10—-71 
follows, and from formula (11) 


dj 
(2) = 0,66. 


Turning now to the experimental data we shall apply the latest 
results of the measurements of Bury and Jenkins, quoted in the Tables 
Annuelles des Constantes et Données Numériques, 1941. The following 
table gives the data for the six lowest concentrations of nitrobenzene 


in benzene: 
A Tia AT 


In the second column are given the values of the lowering of 
the freezing point in an ideal solution, calculated by formula (1), in 
the third column — the measured lowerings, and in the fourth co- 
lamn — the values j calculated by formula (2). The course of j as 


nitrobenzene 


a funetion of m is represented in Fig. 2. It can be seen that the tan- 


: dj 0.66 
gent at the point m=0, calculated theoretically, i. e. rm NR; 
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fits fairly well with the experimental curve. Experimental data for 
concentrations lower than m=0,1 are, however, lacking. 

To verify that the relation between j and yu resulting from the 
integral (16) fits well with the experimental results, the theory has 
been applied to solutions of other polar substances in benzene, viz. 
to solutions of ortho-nitrophenol and para-nitrophenol, as well as to 
ortho-crezol and para-crezol. 


70H is 
1. Ortho-nitrophenol, | p=3,1.10-8, a=—,=0,250.10-21, 
4 kT 


For the minimum possible distance between the dipoles of two coupled 
molecules, the same value is assumed as for nitrobenzene, namely 
ra= 4,6 A, as this distance is determined by the largest group NO, 
common to both these substances. From these data, graphical inte- 
gration gives J=0,50.10—#!, and consequently (3) = 0,13. 
Mj m=0 

In Fig. 3 it is seen that the theoretical tangent fits very wel 

with the experimental curve (these and the further experimental data 


are from the Tables Annuelles des Constantes et Donnée 


1941). s Numeriques, 


NO, 


2, Para-nitrophenol, #=5,02.10-8, a= 0,657.10—1, 
ò 


ro=4,6 A, consequently J=8,0.10-2! and (2) =2,12. It can be 
Mj m=0 i 


~x ES l 
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iseen that the increment, of the dipole moment from 3 0019) og) for 
‘orthonitrophenol, to 5, 02. 10—!8, as for paranitrophenol, gives a 16-fold 


dj 
jincrease in (72) , Whereas the experimentally found increase is 
m=0 


‘somewhat greater. 


OH 
3. Ortho-crezol. an =1,41.10-8, a=0,0518.10-2!. In 
Ne 


j 


© eli Seve elke sey." 
Fig. 4 


this case, the minimum distance between the two coupled molecules 
would be much less than for the foregoing substances. A near enough 


© Q 0A Ole OF 
Fig. 5. 


conformity with experiments can be obtained (cf. Fig. 4), by taking 
ra=1,98 A. Then J=0,61.10-# and (3 ial = 0,16. 


m=0 
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CH, 3 
4. Para-crezol, © p#=1,57.10—-8, a= 0,0642.10—". Assuming 
0 - 


uy) _ 
the same r as for ortho-crezol one finds J =1,63.10—! and ee a = 
M} m=0 


= 0,43; this last result is almost exactly verified by experime 
(Fig. 5). . 
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ON NON-LOCAL QUANTUM ELECTRODYNAMICS 


By Jerzy RAYSKI, Department of Physics, Nicholas Copernicus 
University, Torun. 


(received April 11, 1951)* 


A general rule for constructing observable densities enables the formulation 
of a non-local electrodynamics. The basic equations describing the interaction bet- 
ween the electromagnetic and the electronic fields are integral (or mixed integro- 
differential) equations. The free waves and the non-local analogs of Green’s func- 
tions AT, 424% may be constructed by invoking Born’s principle of reciprocity 
(i. e., a symmetry between position and infinitesimal displacement operators). The 
S-matrix may be constructed by Yang’s method. The commutation relations for 
the (perturbed) fields may be computed from the field equations, provided the homo- 
geneous integral equations possess no solution. The conservation laws are secured 
in spite of the absence of continuity equations. Constant electromagnetic potentials 
are shown to be physically meaningless. The non-local formalism is free from the 


usual convergence difficulties. 


1. Introduction 


In spite of a considerable success of the modern quantum field 
theory, developed by Tomonaga, Schwinger, and others, the funda- 
mental inconsistencies and difficulties have not been satisfactorily 
solved but rather (partly) avoided. Schwinger’s idea of renormalization 
has proved powerful but the renormalization technique is certainly 
mathematically incorrect so long as the expressions to be renormalized 
are infinite or ambiguous. The idea of formal regularization, exten- 
sively discussed by Pauli and Villars, avoids the dangers of ambi- 
guity of the mathematical expressions encountered in the course of 
the computations but cannot be reconciliated with other fundamental 
principles of the local quantum field theory. On the other hand, the 
very attractive idea of realistic regularization, in spite of a conside- 
rable success in removing divergences of the mass-renormalization and 
‘self-stress types, has proved uneffective in the problem of self-charge. 
. In view of this situation it is felt that no half-remedies (more 
or less skillful) are sufficient but we have to look for a radically new 
‘approach to the problems of fields and elementary particles, maybe 


* Revised manuscript received February (17, 1952). 
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as revolutionary as, years ago, was the theory of relativity or quantum 
mechanics. These formalisms have had a remarkable feature in com- 
mon: both consisted in an incorporation of a universal constant (c or h) 
into the theory. The opinion has been often expressed that the reason 
for the inconsistency of the field theory is due to the omission of an- 
. other fundamental constant 4 with the dimension of length which has 
not yet been observed in experiment owing to its extreme smallness. 
Expressed in natural units: e=h/2x=A=1 all physical quantities will 
become dimensionless, which gives us a hint that the incorporation 
on a fundamental length 4 will be the last step completing this period 
of development initiated by the relativity and the quantum theories. 

Born’s and Yukawa’s ideas of non-localizability offer a vast 
generalization of field theory and provide us with a considerable 
amount of hope that in the framework of a non-local theory it will 
be possible to introduce consistently the conception of a fundamental 
length 4 and to solve the chief difficulties inherent in the traditional 
local quantum field theory. 

It must be remembered, however, that the incorporation of uni- 
versal constants into the theory was accompanied by sacrifices of 
some of the customary conceptions and by renunciations of some no- 
tions and postulates which were formerly regarded as obvious and 
indispensable. Just as relativity compelled us to abandon the notion 
of an absolute time and quantum mechanics the old conception of 
state and measurability of canonically conjugated quantities, so the 
introduction of a universal length must violate some of the traditional 
conceptions. The author believes that the incorporation of a universal 
length must be accompanied by the renunciation of differential equa- 
tions (as means of description of field interactions). The differential 
equations have to be replaced by some integral equations which enable 
a relativistic kind of an action at a distance. In consequence it will 
be impossible to retain the old conceptions of densities obeing con- 

tinuity equations as it is obviously impossible to trace in detail the 
currents and the energy and momentum flow inside the domain of 
the order of A, that is, inside the elementary particles. This resigna- 
tion does not constitute a decisive difficulty since it is possible to 
guarantee integral conservation laws in spite of the absence of con- 
tinuity equations. = 
In a previous paper (Rayski 1951), denoted heraafter as (I), the 
ideas of non-localizability and reciprocity were applied to the (rather 
academic example of) two coupled scalar fields. The present paper is 
chiefly concerned with non-local electrodynamics, but contains also 
some improvements of a general nature. 
— d 
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2. Observables 


The idea of non-localizability consists in the assumption that 
field variables depend not only upon the operator of position but also 
upon the displacement operator 


p=Yy(X,, d,), (1) 
where 

EF d, |= ôw. (2) 
It will be convenient to introduce the following (dimensionsless) com- 
binations of a, and d, 


EE paas 1 = gz 
at= z (2 ‘at + ZA du), a, = (2 *aty—> dy] (3) 


[at a, ]= ibu, (4) 


where A is a constant with the dimension of length. 
The transformation 


>ii, Ady ida, . (B) 
yields 
SC SS (5’) 


and leaves the commutation relations (2) and (4) unchanged. 
Introducing a representation in which «, are represented by 
diagonal matrices the field variables y assume the form of matrices 
<x'/p/a’’>, where Su, #, denote (continuous) eigenvalues of the oper- 
ators Xy. In order to avoid confusion with representations connected 
with field quantization (which play a somewhat different role), it will 
be convenient to introduce a tensor-like notation for such matrices. 
We agree to label y with lower and y* with upper indices: pwe, p**’*” 
The product of a field variable and its conjugate is the matrix product 


Ih yl ot 
on = fy" AL Puig) (6) 


where da denotes, in short, the integration over the fourdimensional 
continuum. The formula (6) suggests that observable quantities are 
attached to matrices with an upper and a lower index. Therefore, we 
agree to denote matrices describing real (i. e. observable) fields by 
means of one upper and one lower index. For example, we may os 
note the matrix describing the electromagnetic potentials by (Ajy)y , 
and the condition of reality becomes the condition of hermiticity 


(Ane (Aya. (7) 
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Possible representations for x, and d, are 
At f] 1 
(Suy = 0 O(a’ —2""), (dy) =— sr ôe — e), (8) 
OL 

where 6(a) is the fourdimensional Dirac delta function. By introducing 
Yukawa’s variables 
Lut bu 


9 ? 


p 


A= Tu Lu Las (9) 
the field quantities may be regarded as functions of these two kinds 
of variables: y=y(X, r), where X, are connected with the „centre of 
the particle” while r„ are connected with the internal structure of the 
particle. 

In order to introduce an interpretation of the non-local theory 
we define a hermitian operator I(x) whose representation is the four- 
dimensional Dirac delta-function 


Tei =o" —a), (10) 
With the aid of this operator we establish the following interpreta- 
tion 1: Let æ be a hermitian operator constructed of field quantities 
of a (single) field, then the corresponding observable density a(a) is 
given by 

alx)=Tra I(x). (11) 
The definition (11) establishes, at the same time, a correspondence 
with the local field theory, where a(x) is constructed from the local 
field quantities essentially in the same way as @ is constructed from 
the corresponding non-local field quantities. If, for example, a(x) means 
the charge and current density j,(x”) for a spinor field, then in the non- 
local theory we have 


gs = ie fp x" yl! RE yy (12) 
and 
jv j= LTD I(x KAE dis ALT T(x))*), dar’ 


=f [aX dri X, 1)(X—2)= får j (æ, r) 


where the point «, is a substitute for Yukawa’s variable X, given 
by (9). More generally, if æ is of the form 


(12") 


Buy... =P Duv.. P, (13) 
t The same interpretation has been given independently by Meller and Bloch 
(Bloch 1950). 
2 We shall often write xu instead of Xy. 
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where Oy»... is an operator constructed of Dirac matrices, we find 


p =f f fr o dx" pyr uda" (I(@) xr ; dæ’ 
x£ , 7 vit 
ire =| d(x anae a gee! m a) Ope 


‘oad nae 


(13°) 
2 


=f far dry p[2—5, T Jow ple- 5, r), 


where 7’ =a'—a'"", r’=a"'—a'". Thus, the correspondence rule for bi- 
linear densities written in terms of position and structure variables is 


PE) Ow.. Y w> f far dv y(o—F, r a v(o— D zor’) (14) 


3. General form of field equations 


Our next task is to generalize the equations of quantum electro- 
dynamics 
O Ault) =— juli) (15) 


(14 ge +) ple) = iey“) Ana). (16) 


Equation (15) establishes a connection between the observables A, 
and j, so that the interpretation rule (11) immediately yields a plau- 
sible generalization 


oTr A, Iæ) =—Tr ju I0). (17) 


This equation may be integrated by means of retarded potentials 
which yields 


A,(v)=TrA,I(«) =Tr A; I(x) + [de' D™e—a') Trin dpa (18) 


A similar solution may be written in terms of advanced potentials. 
The solution (18) may be substituted into (16). We see that only the 
observable ,,density” ea Gees I(x) but not the non-local matrix 
(Ay et appears explicitly in the field equations, whence it is to be 


inferred that, for real fields, the non-localizability is only apparent. 
Thus, we may simply ite always A,(#) instead of ZrA,I(x) and 
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treat the electromagnetic potentials as localizable. Equation (17) and 
its solution (18) are then simplified to 


DA,(#)=—Trj,L(x) (17°) 
A,(v)=Ap(2)+ [da'D™e—a')Tr ja Ila’). (18’) 


The problem of a generalization of the Dirac equation (16) is 
more difficult since we have to do here with essentially non-local quan- 
tities y. Let us integrate (16) by means of retarded (or advanced) 
potentials: 


vio)= ya) +5 f da’ Swa yiye Ayla). (19) 


In this way the differential equation has been replaced by an integral 
equation which is more suitable for generalizations introducing a non- 
local point of view. We assume that non-local counterparts of the S 
and 8 functions exist, and introduce formally their combinations 


set—§—19, swv—S+19. (20) 


We postpone to later sections an explicit construction of these func- 
tions and limit ourselves, for the present, to the remark that these 
quantities are operators whose representations must be of the form 

ATHENA 5 ` 5 
Syvx with two lower and two upper indices (since they are closely 
connected with the diadic product of two field quantities y y, and 
p=’). These matrices may be regarded also as functions of Yukawa’s 

variables 

xe’ + a" g” + glV 

2 


ye) Ae , 
SAI O GET; 7 


SKO al —gIV— ¢. (21) 


From the requirement of invariance with respect to space-time trans- 
lations it follows that these functions must depend only upon the 


difference of the position variables X®—X® =X. Thus, we may write 
SSC r 2) SS Xr), (22) 


Using the correspondence rule (11), a plausible generalization 
of (19) is 


pares faw PS 744 Ta) Ayla’) (23) 


This integral equation may be written more explicitly in the matrix 
form 


(23’) 
in € , 11 i 
Pare Vx TF af da i J ip (Stet yt da y rv da’ (I(x) E}, dx” Ala) 
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or, going over to Yukawa’s variables, 


pla, URAR r) 


analogous formulae hold for the advanced potentials. 

Reassuming, it may be stated that (i) the observable densities 
are given by means of the general rule (11), (ii) the basic equations 
of the formalism are the integral equations (18’), (23), and similar 
equations written with the aid of the advanced functions 8%4, Dadv, 
The electromagnetic potentials may be treated as localizable without 
loss of generality. The quantities S'@4”) are given by (20), where S 
and § are operators which, in terms of Yukawa’s variables, may be 
written in the form (22). The formalism will be completed by con- 
structing free ingoing (or outgoing) waves and the functions S and 8. 


4. The reciprocal S function and the interaction-free field 


By definig a suitable S8(#,r,s) function the commutation rela- 
tions between the free ingoing (or outgoing) wave functions and the 
free field equations will be fixed. 

The principle of non-localizability which was exploited in Sec- 
tions 1—3 cannot determine the form of the S function unambiguously. 
To construct an S function a further principle is needed. To this end 
we invoke the principle of reciprocity. In conformity with this prin- 
ciple a A(a#,7r,s) function may be introduced in the following way 


Alæ, r,s) sau f ipera ,7,8), (24) 


where the Fourier transform is 
A(p,7,8)=N(p) o(p,7) A(p) e(p, 8) N*(p) (25) 


and A(p) is the Fourier transform of the usual 4(x) function well 
known from the local field theory: 


A(p)= — 22 106(pi, + m?) sgn Do, (26) 


while o(p,r) is short S the product of two (one dimensional) delta 
functions 
oe (Lat(r, p)P) (Laz (7, p) P). (27) 
The quantities . } 
a(r, p) =A TuE App (27°) 
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are closely connected with the operators (3). N(p) is a normalizing 
factor so that 


N(p) f olp,r)ar=1, N*(p)=N(p). (28) 


From (25), (26), and (28) we have a simple connection with the tra- 
ditional A function 


f [eras A(a,r,8)=A(a). (29) 


From (26) it is seen that the new A-function satisfies the Schrédinger— 
Gordon—Klein equation in terms of the position variables 


(> —m!)A(a,r,8)= 0. (30) 
u 


Owing to the existence of other 6-functions in (27) some further rela- 
tions, namely 


i 3 \ J2 3 
Gat =x) A(v,r,s)=a#(s, : =) A(u,7,s)=0, (31+) 
) 


are satisfied. Taking account of (31 
also in the form 


equation (30) may be written 


a LEIN ape gh ac 
ag (r, a 5a} Oy (r, 7 5g) 40 s)= 2m? A(x,7,8). (30’) 
The same equation holds with r, replaced by s,. 


We notice the following symmetry properties of the generalized 
A-function 


A(— Lu? yp )8u) =—A(Lp, 1p, Sp) (32°) 

A(Lu,—T us Su) =A( Ly Tu) —Su)=4(Lu, Tu; 8p) (32°) 
-=> => > >> > 

A(x, T, 8—20, — To, — 80) =—A(2, r, 8, Ly jo, Sq). (32) 


Thus, 4 is skew-symmetric with respect to the transformation 2,>—# 
symmetric with respect to the transformations r,+>—r, or MEE 
and skew-symmetric under the reversal of the time axis pices 
To>—fo, So™>— So. It should be noticed, however, that the golieralifed 
A-function is not skew-symmetric in 2) alone and, therefore, does not 
vanish identically outside the light cone of the position vaia £ 
This new property is intimately connected with the non-local shee 
racter of the field. 

With the aid of the above generalization of the A-function we 
define the generalized S-function by means of 


: l 
Sagl, r,s) =h z m) Ale, r,s) (33) 
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‘and postulate the following commutation relations for the free (in- 
i going or outgoing) waves 


{pigeon peo 


1 
Lr), Y x $)}= 7 Sap(e—2"', 1,8) (34) 
while all other anticommutators vanish. The same commutation rela- 
tions hold for half the sum of the ingoing and outgoing waves 


serene ne iver : 
{pO(x,7), y(x S)J=— Supl — E 8), ye Yt yout) (34°) 


The anticommutation relations (34) or (34’) are satisfied by free wave 
functions of the following form 


yO(a,n) = Ae salts YOP) N(p) o(p,"), (35) 


where y(p) are spinors well-known from the local theory of fermions 
in vacuo, while w#(p) are the usual annihilation (and creation) opera- 
fors for negatons(and positons). Here again o(p,r) and N(p) are given 
by (27) and (28). The free waves (35) satisfy the same equations as 
the S-function, that is, the Dirac equation 


| o,r) 0 (36) 
and the equations 
10 1 ə 
Ga @ m z) Oe ("5 =| yO(a,r) = 0. (37+) 


As a consequence of the Dirac equation the non-local spinor wave 
tunctions satisfy the Schrédinger-Gordon equation which may be writ- 
ten also in the form 


eee fae 
ap [3 an G 7 A pO, r)= Pm, r) (38) 


owing to’ (STE). 

The commutation relations and the field equations for inter- 
action-free spinor fields have been constructed thus far in the w-re- 
presentation, but we may go over also to the operator form where 


equations (37+) and (38) become 
[at, [ar, p]]=0 (374), [a7, (a7, pOl]=0 (371) 
fat, (az, pO = my, (38’) 


where the operators ax are given by (3). The same relations hold for 
the ingoing and outgoing waves separately. The equivalence of (37') 
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with (37+) and the equivalence of (38’) with (38) follows from (3), 
(Z7’'), and (8). ; 

The equations (37’) are obviously reciprocal since they inter- 
change under the transformation (5°). It may be asked whether equa- 
tion (38’) is reciprocal or not? From (4) it follows that 

[at, (a>, f=le,, La, fI (39) 
for an arbitrary f. In this way we see that the transformation (5) 
changes the sign of the left-hand side of (38). Thus, (38) is not in- 
variant under the transformation (5) but it is invariant under a trans- 
formation 


Aa, —>tAdy, Ady >is ay, Am—>idm (5) 


whose meaning, however, is not altogether clear. 

Concluding it may be stated that free waves obey the completely 
reciprocal equations (37+) and the Dirac equation. The Dirac equation 
is not reciprocal, but equation (38) which is a consequence of the 
Dirac equation is reciprocal in the generalized sense: it is invariant 
under the transformation (5”). 

We cannot expect a reciprocal generalization of the Dirac equa- 
tion itself since the reciprocity establishes a symmetry between po- 
sition and displacement but not between position and the rotational 
degrees of freedom connected with Dirac’s matrices y,. 

A further remark may be added to clarify the question of uni- 
queness of the solution. By postulating the form of the commutation 
relations we find the unique solution (35) for the interaction free field 
which obeys three equations: Dirae’s equation and equations (37+) 
(while (38) is only a consequence of them). It would not be so if we 
proceeded in the oposite direction, that is, if we postulated the Dirac 
equation and the equations (37+) and tried to infer the form of the 
solutions and the commutation relations. The three above equations 
are not sufficient to determine completely the solution. Instead of (35) 
we would get a solution containing an arbitrary (though regular) factor 
dependent on the internal structure variables r,. This factor intro- 
duces a possibility of particles with higher spin values discussed by 
Yukawa (1950) and Fierz (1950). We purposely drove in another di- 
rection to get a theory of particles with spin 1/, only. The lack of uni- 
queness of the solution of equations (37+) and Dirae’s equation may 
be seen from the fact that, in the framework of the local theory, we 
have five equations instead of three. These are the Dirac equation 
and the four conditions of localizability [Xu p]= 0 for w=1,...,4 
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5. Further discussion of observables 


The two delta functions which appear in the solution (35) may 
ibe replaced by means of (27) by two other delta functions 


o(@,7r)=10(A ee: aye Du) Tn Pu), (40) 


‘which shows that the relations (37’) may be replaced by the equivalent 
irelations 


A Hu, [Eu, y]]= Mi du, [d., vii, [Xu [dy, y]]=0 (37) 


iused in (I). 
Owing to the existence of 6(r,p,) we have, from (35), the fol- 
lowing identity 
p(“&—kr,r)=y(x,r), (41) 


where & is an arbitrary constant. An analogous identity holds for the 
S-function 
S(a—kr—a’+ k’s,r,s) =S8(x—2’,r,8). (41’) 


This enables us to write rule (14) for constructing densities in a sim- 
pler form. Denoting in (14) r’ by r and r” by s and putting in (41) 
k=+4 the rule for the transcription of densities into the non-local form 
becomes 


(2) Ou: ple) | | drdsiply,s) Our..0(Y,"), (42) 


where 
y=æ— (rF s). (43) 


This rule may be stated as follows: in order to obtain a density (being 
a bilinear form in y and wp) at the point x replace the local quantities 
p(x), p(x) by y(x,r), p(a,s), shift the point x to y given by (43), and then 
integrate over d'rdts. 

It is easily seen that for an interaction-free field the new den- 
sities satisfy the continuity equations and yield the same eigenvalues 
for the integrated observables as in the local theory. Let us consider, 
for example, the charge and current density 


jula)=te f fardsyly,s)y*yly,”). (44) 


The non-local field obeys the Dirac equation which may be written 
also in the form 


(rg +m) vy) =0 (45 i 
“M 


ea 7h mye) Uk ae) 
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by multiplying (45’) by y(y,s) to the left and (45’’) by y(y,7) to the 
right and substracting one result from the other, we get 


vA) ye 
Bn, P UAE 0, (46) 


whence the continuity equation for the charge and current four-vector 
given by (44) follows immediately. 

The total charge is obtained from (44) by integrating jọ over 
the threedimensional space 


[@xj=ef@af fardsply,s)vy,) 
=o X [Pa f fdrdsctoovN*(q) o(9,8)0*(qy*(q)v(p) Up) ops) Np) (47) 


Pq 
=e X f [ards N*(p) o(p,8)u*(p)ee(p) o(p,r) N(p)= e Ju” (pulp), 
P p 


where use has been made of (35) and (28). The last expression is iden- 
tical with that in the local theory and, of course, may be charge-sym- 
metrized. 

In the same way the continuity equations may be inferred for 
the generalized energy — momentum — stress tensor and it may be 
shown that the operators of energy and momentum are identical with — 
those in the local theory. Thus, in spite of,a slight difference in den- 
sities, this non-local theory of electrons in vacuo yields identical results 
with the local theory. 

It may be shown that, in the limit A=0, the interaction-free 
waves, the commutation relations, and the densities go over into the 
local ones. For free particles the four-vector p he is restricted by oo — mè. 


Hence we may consider any plane wave picked out from the solution 
. . = l 
(35) in its system of rest p=0, pe=+m. The two delta functions in 


(27) or (40) yield r= 0 and amt, wherefrom we see that the values 
of the structure variables r, for which y(#,r) is different from zero 
tend (separately) to zero for A+0. This is true in any frame of refe- 
rence. Hence, for 4+0 the vectors r, (and s,) for which y(a,r)=:0 
vanish so that the variables y, may be replaced by a, in the densi- 
ties (42). Then, taking account of (28) and (35) it is seen that the 
densities reduce to the local ones. In this way the correspondence 
with the local theory is fully guaranteed. 


6. The S-function 


Let us return to the problem of the interaction of electrons with 
the electromagnetic field given by means of the integral equations (18’) 
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and (23) or analogous equations written in terms of advanced poten- 
tials. The ingoing (and outgoing) waves appearing in (23) are already 
known since they are the interaction free waves discussed in Section 4. 
Sr and S are defined by (20) where S is given by (33). We do not 
know yet the generalized Green’s function §. Since the spinor field 
and its interaction with the electromagnetic field are essentially non- 
local, we cannot expect the generalized § function to obey the tra- 


- ditional equation (rage +m) = —ò(x). For a non-local theory the 
ub 


right-hand side (as well as the left-hand side) must depend upon the 
structure variables r, and s, so that it cannot possibly be a delta fun- 
ction of the position variables æ„ alone (which was a characteristic 
feature of the local theory reducing the integral equations to differen- 
tial equations). The non-local theory is expected to involve integral 
equations that are not reducible to pure differential equations. 

A plausible form of an S(#,7,s) function may be constructed 
by invoking again the principle of reciprocity. We demand that S 
satisfies the reciprocal equations (37+). Such an S function may be 
constructed on closely analogous lines as the S function 


s 2 
Sao mS) ( naa Alans) (48) 
where the Fourier transform of the A function is j 
A(p,r,8)=N(p) o(p,r) A(p) o(p, 8) N*(p), (49) 


where A(p) is the Fourier transform of the usual 4 function well- 
known from the formalism of Schwinger: 


A(p)=P.v. (50) 


while o(p,r) and o(p,s) are given by (27). The o(p,r) and o(p,s) factors 
guarantee the validity of equations (37+) also for the 4 function which, 
in turn, ensures that not only the unperturbed functions y'*, y™, or y® 
but also the perturbed y given by the integral equation (23) satisfies 
the reciprocal equations of the form (37+). It goes without saying 
that the perturbed wave function y does not satisfy the relation (38). 
Instead of (38) (which is valid only for unperturbed fields) the per- 
turbed yp satisfies the integral equation (23). 

l A considerable complication arises from the fact that (contrary 
to the case of the A function) the variables p, appearing in a A(p,r,s) 
are not restricted to p= ne. For space-like p "A the integral 
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is ambiguous and, therefore, the definition (28) of the N(p) needs 
a refinement. In en 10 we shall give a more complete definition 
of the normalizing factor N(p) which reduces to (28) for time- like py. 


7. Formal solution of integral equations 
Heisenberg’s S-matrix 


From (11) and (44) the solution (18’) may be written in the form 


AN ae gee i AP'D*(@—e' pus), yry (52 


where 
r+ e’ 


dP'=dx'dr'ds' = dy'dr'ds', y'=% — — 


(53) 


Here the charge and current have been written, in a symmetrized 
form. In the same way the integral equation (23’’), describing the 
influence of the electromagnetic field upon the spinor field, may be 
written in the form 


in e gre , 1 DE-A? ’ 
v= +s f AP'S y yrs O) Ada}, (54) 


where the free variable y has been substituted for #, and æ'—4s' or 
æ'—+4r' could be replaced by y’ owing to the identities (41) or (41’). 
Moreover, the product yA, has been written in a symmetrized form. 
This symmetrization is necessary since the perturbed non-local quanti- 
ties y do not commute with the perturbed A,. Except for a symmetri- 
zation, the equations (52) and (54) are equivalent to the original equa- 
tions (18’) and (23). Quite analogous equations hold for advanced 
solutions and for outgoing waves. 

The S-matrix may be computed from the integral equations by 
means of standard methods developed by Yang and Feldman (1950). 
The formulae for the computation of the S-matrix elements are 


Ay(a)= AQ a) +> SAP Dao BYS h (55) 
gee dP'S(y—y',7,8') {ypy r), Aula"), (56) 


ENO s=- f aP'D(x—2') {S [p(y’,s’), y pyr] (57) 


WPu), S=- f aP'S(y —y',r,s') {S, {y"vly',7'), Aula’)}}, (58) 


where 40 and y® mean half the sum of the in- and outgoing waves 
and obey the commutation relations (34’). For the derivation of equa- 


oe 
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tions (57) and (58) see, for example, the previous paper (I). From 
(34) and (55)—(58) the S-matrix may be computed by iteration but 
the result may be guessed immediately: since solving equations by 
means of iteration is a mechanical procedure, the new S-matrix ele- 
ments will differ from the traditional ones by exactly the same sub- 
stitutions which convert the corresponding formulae of Yang and 
Feldman into our formulae (55)—(58). Thus, in order to get the S-ma- 
trix elements for the non-local theory we have to substitute in the 
traditional S-matrix elements (derivable directly from Yang equations) 
every p(x) (and p(a)) by py, r) (and p(y, s s)), every S(x — g") 
(or A = by S(y’—y",r’,s’) (or by S(y’—y”,7’,s’’)) while the 
integrations dæ'dæ”... have to be replaced by threefold integrations 
aP’dP’... given by (53)? 

Thus, the existence of an S-matrix is (formally) proved by means 
of its explicit construction. The (formal) existence of the S-matrix 
constitutes, at the same time, a proof of the compatibility of our 
integral field equations among themselves and their compatibility with 
the commutation relations (34). The real proof of the existence of an 
S-matrix will consist in showing that its elements are free from the 
usual divergences. 

If only the inhomogeneous integral equations (55) and (56) pos- 
sess a solution then we may find the perturbed field y from (55) and (56) 
by iteration. The anticommutation relations between the perturbed 
non-local y, y functions follow automatically from the integral field 
equations and from the commutation relations between the unper- 
turbed functions y® and y®. These commutation relations may be 
computed, together with y and y, by iteration and expressed in a power 
series expansion of the coupling constant e. Thus, we possess a direct 
method of quantization of the field equations completely avoiding the 
detour through the canonical formalism. 

The situation will be considerably complicated if, under circum- 
stances, the homogeneous integral equations bd from (55) and 
(56) a dropping the free waves y® and A®) possess solutions. The 
author believes that such extra solutions exist and describe stationary 
states that cannot be obtained by means of the customary pertur- 
bation treatment. Special methods have to be developed to solve the 
homogeneous equations. Probably it will be possible to orthonormalize 
these solutions and to attach to them (independent) creation and_an- 
nihilation operators ing order to complete the second quantization. 


_ 8 When solving the equations care must be taken of the proper order of fac- 
tors since Ay are not taken at exactly the same points as yp. 
Acta Physica Polonica 9 
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This problem has not yet been solved even in the framework of local 
electrodynamics. 

In general a solution will be extremely complicated but, from 
physical arguments, we may expect that the homogeneous equations 
will possess no solutions unless a third field (e. g., a proton field) is 
involved, or unless A, is a suitable external potential. The reason 
is that there cannot exist true stationary states for electrons alone 
(positons bound to negatons must annihilate sooner or later). 


8. The auxiliary condition 


The problem of the auxiliary condition is very simple in the 
framework of a formalism expressed in terms of integral equations. 
The perturbed potentials A, are given explicitly by means of (55) 
and may be eliminated entirely by introducing (55) into (56) and (58). 
In this way the field equations (56) and the expressions (57), (58) de- 
fining the S-matrix are expressed in terms of the unperturbed poten- 
tials AR only, while (55) may be regarded as a definition (of a second- 
ary quantity) without influence upon the fundamental equations. 
Therefore, we may limit the auxiliary condition to unperturbed po- 
tentials, 

Fay AM te) |> =0 (59) 
being a restriction upon the state vector |X. This restriction reduces 
to two well known initial conditions as usually, owing to the fact that 
A® satisfies the wave equation 


oA —0. | (60) 


9. The problem of gauge invariance © 


The non-local theory is clearly invariant with respect to gauge 
transformations of the first kind 


pla, r) > (a, re“, p(w, 8) > pla, 8) e (61) 


(with a constant a) since the observables of the non-local Spinor wave 
functions are always bilinear in y and its conjugate y. On the other 
_ hand, the above formalism is not invariant with respect to gauge 

transformations of the second kind. The author does not share the 


< Ti 
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‘opinion of (a great majority of) physicists who consider the gauge 
iinvariance of the second kind for an established and indisputable 
property of electrodynamics. Several doubts may be raised against 
the postulate of gauge invariance of the second kind. First of all, we 
‘cannot be sure that Maxwell’s equations are exact. The very fact 
that it is impossible to do without introducing potentials explicitly 
(they appear in the Lagrangian and the Hamiltonian) makes it plau- 
sible that in a future development of electrodynamics not only the 
field strengths but also the potentials A, themselves or their deriva- 
tives may acquire a certain physical meaning. Moreover, it is not sure 
whether the rest mass of the photon is exactly zero; it may well as- 
sume a finite though extremly small value. But even if the photon 
rest-mass is exactly zero, the gauge invariance of the second kind is 
not a necessary but only a sufficient condition for its vanishing. A value 
zero for the photon rest-mass may be secured by other devices, e. g. 
by a realistic compensation or by a formal renormalization, and not 
necessarily by the gauge invariance. It may be stated quite generally 
that the problem of the particle rest-masses lies beyond the scope of 
the present-day theories. The photon rest-mass seemed to be the only 
exception in this respect but such exceptions are suspicious. We prefer 
therefore to treat the rest masses of all kinds of particles including 
photons on the same footing, i.e. adjust their values by means of 
renormalization to agree with experiment. 

At any rate, we have to guarantee a covariant form of the equa- 
tions with respect to phase transformations 


p(x, r) > p(x, r) ea (62) 


with an arbitrary A(w) since the observables of the non-local theory 
(42) are bilinear in the wave functions y and y taken at the same 
point y given by (43), so that phase transformations of the type (62) 
must not influence the physical content of the theory. 

Let us introduce generalized equations (37+) 


12 eit ee: 
(ra + = 3p, —Anler) yala, =O, (63+) 
1 ə 2 
(nF dar, + Ay(2) pa(@,7r) = 0, (635) 
where | 
Ail) = xe A) (64) 


Q* 
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(We are using now units c=h/2x=A=1). The equations (63+) are 
clearly equivalent to (37+) since their solution y, differs from that 
of the equations (37+) merely by a phase factor 


pa(@,r)=Yy(£,1) Pagid! (65) 


Of course, the change of phase (65) affects at the same time the Dirac 
equation in vacuo and the commutation relations: 


ə . ] 0 PEET f (0) na zn (0) m! N Ges ) (66) 
Ve a to) pesca Wa (x, r) po (x S) =S (0,0 7,8), 


Wy 
where uns 
Sala, a’, 7,8) == eAMS(~—a’, 1, 8) eA 
ene gn (67) 
= |y“ or ae: —m|Ap(x,x',7, 8) 
with nie 
Aale, e, 1,8) = e4@) A(a—y’ r,s) e-4@). (68) 


lt may be easily shown that the phase transformed wave function pa 
obeys the following integral equation 


e Dei , / Donat 1 
vas = NH 5: f AP Bayr, Spay) Ade} (69) 


if y(y,r) obeys equation (56). Equations (63+), (66), and (69), together 
with (55) (which obviously may be also written in terms of y4), are 
equivalent to the analogous equations with A4=0 but their advantage 
consists in the explicit covariance with respect to phase transformations 


par=pae#), A'(2)=A(aw) —a(a). (70) 
The phase covariant transcription does not introduce any new 
physical content, and is an analog of the transcription of the equa- 


tions of the special theory of relativity from their standard form (in 
Lorentz coordinates) into a covariant form valid in general curvilinear 


coordinates. The fictitious potentials Ay) = Aa) play a role similar 
tt 


to the fictitious gravitational potentials provided by general coordi- 
nates in a flat space. Those fictitious gravitational potentials may be 
removed by means of a suitable coordinate transformation, while the 
fictitious potentials A, may be removed by a suitable phase trans- 
formation. 
Once the covariance of the formalism with respect to phase trans- 
formations is established, we do not need to introduce the true electro- 
magnetic potentials A, at the same places which are occupied by the 
fictitious potentials Au, provided we are ready to ascribe to the Ans 
or to their derivatives a certain physical reality. 
Sa . 
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We now show that even in the framework of our non-gauge- 
Invariant theory constant electromagnetic potentials A „=Q, are physi- 
cally meaningless since the interaction with such potentials produces 
merely a change of phase. Let us apply the Dirac operator to (56) 
and replace A, by the constant potentials Cu, 


ye bm p= iE; —me) fara Fins fe 
OY ONY ot a dy? m A(y—y 7,5 ) {y yy Y ), Cu} 


= 1 (2 pikty—y) 
— 7 Cua (z =m?) fap fonw elkr) Ta elk 8 * l k)yly r) 
u 


e Ji 
a i I dk X ety eow N(k) o(k,r)olk,s’)N*(k) 
Pp 
ulp\ylp) olp, N (p) 


=F, | faras far Seto a(p—b) N(R) ofr) olkys") Ck) 
p 
ulp)ylp) olp, r )N(p) 


e à i 
=— Fyn | faras S Np olp, op, )N*prulp)up) elp, Nip) 
p 7 


=—F "Cu X ePveelp)y(p) elp, r) Nip) 
P 


e 
ars Ths Ze CuvlY, Hye 


Replacing y by x, we have 


[r ( —ie Ca == m) ylæ,r)=0 


IL u 


which shows that 
yle, 7) = p(x, 7) ewe, (71) 


Thus, constant potentials may be absorbed by the interaction-free 
waves y® by means of a change of phase. The result (71) shows that, 
in spite of the absence of full gauge invariance of the second kind, 
constant electromagnetic potentials are physically meaningless even 
in non-local electrodynamics. The invariance with respect to changes 
of potentials by constant values assumes an intermediate position 
between gauge invariance of the first and that of the second kind. 


10. Discussion of the results 


In Section 5 we have shown that the observable densities ob- 
tained with the aid of rule (42) satisfy the continuity equations in 
the case of an interaction-free field. This is no more true for the per- 
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turbed fields whose interaction is essentially non-local. In spite of the 
absence of continuity equations, the integral conservation laws are 
still satisfied. The conservation of charge follows from the invariance 
of the non-local formalism with respect to gauge transformations of 
the first kind while the conservation of momentum and energy for 
collision processes follows from the invariance under translations. But 
this may be seen more directly from the form of the new S-matrix 
which yields, in the same way as in the case of the traditional S-matrix, 
non-vanishing probability amplitudes only for processes for which the 
total energy and momentum are conserved. A closer examination of 
the new S-matrix shows that the new S-matrix elements differ from 
the traditional ones merely by a set of factors of the type 


R(p,q) =N(p)N" (p,q), (72) 
where 


N-\(p,q) = amen e? o(p,r (73) 


and Pu and qu denote the energy and momentum of electrons for real 
or virtual states. In modern terminology p, and g, denote either wave 
vectors connected with the plane waves y® (p), or arguments of the 
Fourier analysed functions A(p,r,s), A(p,7,s) and A®(p,r,s). The com- 


putation of (73) yields the following result: N(p,q) is a function of 
the argument i 


At 
= 7 (Pea) —p -gj=—* “(Pade “Pv du) (74) 
of the form 
a: 
ANES © for a>0, peso (75’), 
or 
1 =a : 
N(p,q) =, ah —a@ Bc a<0. O 


The normalizing factor N(p) given by (28) is identical with N(p,p) 
and may be regarded as a limit 


N(p) =lim N (p,q). (76) 


‘This limiting process is, sana TA EAA If we perform the tran- 
sition to the limit from the domain a>0 we get, from (75’), 
Np) =lim N(p,g)=+ gaa, PS0 Res ce 


but if we let gp so that a<0 then we get from (75%) 
NAN Deres N(p,q) =9. rh 


a>—0 
4 
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[In view of this ambiguity it seems natural to define N (p) appearing 
in (72) either by means of (77’) or (77) according to whether p, and gu 
jappearing in (73) yield a>0 or a<0. With the aid of this definition 
tthe factor (72) becomes 


sin Ya for a>0 
el, 0) =) as ae (78) 
0 yy a 0 


It is possible to introduce such discontinuous factors (78) into our 
formalism from the very beginning, namely by treating N(p) and N*(p) 
as operators with the property 


1 for a0 
N*(p)N(q) =} (82 a4)? (79) 
0 1» ic) 


It should be remembered that this rule for multiplication must not 
be applied to two N-factors appearing in the same A function. The 


two factors N(p), N*(p) appearing in the same A function form a diadic 
product. | 

The factors (78) strongly diminish the interaction for a—oo. 
Owing to these factors the theory is free from the usual convergence 
difficulties. - 

The case a>0 occurs always when at least one of the four-vectors 
Pu. qv is not space-like. In this case the convergence factor tends to 
unity for 40 as it should do in a theory obeying the principle of 
correspondence. The case a<0 is rather exceptional and may occur 
only in higher orders of approximation than e? where two A functions 
are immediate neighbours in an S-matrix element: 


Aa —2"") Ao!" a")... (80) 


The higher order corrections to scattering processes will be sen- 
sitive to the discontinuous behaviour of factors of the type (78) which 
may provide a possibility of an experimental check of this formalism. 
At any rate, all the effects involving the case a<0 and violating the 
correspondence with local electrodynamics, lie in the relativistic region 
(since small perturbations cannot change the time-like character of 
the four-vector of the energy and momentum of the virtual particles). 
Thus, the correspondence with local electrodynamics is still valid in 
the non-relativistic region. 
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41. An alternative definition of the A function 


We have seen that it is possible to formulate a non-local electro- 
dynamics in which the (perturbed) spinor function y satisfies an in- 
tegral equation (54) or (56) and moreover the two conditions of re- 
ciprocity (37+). Therefore, this formalism may be considered as satis- 
fying the principle of reciprocity. On the other hand, the principle 
of correspondence is limited on account of the o ate vanishing 
of interaction for (virtual) transitions in the region a< 0 with a given 
by (74). It may be interesting to notice that it is also possible to for- 
mulate a non-local theory in full agreement with the correspondence 
principle. To this end we introduce an alternative definition of the 
A(æ,r,s) function, namely 

A(au,r, 8) = —4n(x) A(x,r,8), (81) 
where A(av,7,s) is given by (24)—(28) and n(x) is an invariant gene- 
ralization of the function sgnt: 

1 for # <0, %>0 
ena, OF te pea) (82) 
—1 » @ <0, %<0. 


With the aid of this alternative form of the A function we may 
formulate the integral equations (55), (56) and the equations (57), (58) 
definiting the S-matrix. Since the convergence factors (72) are con- 
nected now with py, or qu appearing as arguments in the Fourier trans- 
forms A(p,7,s), 4% (p,7,s), or yO(p,r) (but not directly as arguments 
of the Fourier transform of a A function) they are now restricted to 
at Se so that the case a<0 never occurs. The correspondence 
is assured since the convergence factors given by the first formula (78) 
tend to unity for A—0. 

On the other hand, this alternative 4-function does not satisfy 
the reciprocal equations (37+) so that the perturbed y given by means 
of (56) does not satisfy (37+) either. In this case the correspondence 
is assured but the principle of reciprocity is restricted to the case of- 


the non-local field in vacuo. It seems impossible to satisfy both prin- 
ciples simultaneously. 
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This paper deals with the theory of the interference effect of the enforced elec- 
tric dipole radiation and the spontaneous electric quadrupole radiation. The case 
is investigated when the magnetic field is parallel to the electric field, the latter 
enforcing the dipole radiation. The result of the calculation is that the interference 
effect can be expected to be observable in transversal observation. Placing a nicol 
askew, we can obtain such a position that one of the AM= +1 components is ex- 
tinguished. Then, the intensity of the corresponding AM=+1 component gets 
its maximum. Approximative formulas for the intensity of the Zeeman components 
are given in the case of Russel-Saunders coupling. Their application to the first 
line of the 3D;/.,3/. > 2S1;, series in the potassium spectrum leads to the conclusion 
that in this case one must apply an electric field of about 10 kV/cm to observe the 
interference effect. 


Introduction 


In a report on multipole radiation by A. Rubinowicz (1949) an 
interference effect in the Zeeman components between the electric 
dipole radiation and the electric quadrupole radiation has been assu- 
med. If one wants to break the Laporte rule, which forbids the simul- 
taneous appearence of spontaneous electric dipole and quadrupole 
transitions, one must apply an external electric field. This field causes 
that besides the already present spontaneous electric quadrupole ra- 
diation there appears also an enforced electric dipole radiation. To 
make such an effect appear clearly, we can take advantage of the 
enforced character of the dipole transitions regulating its strength by 
changing the enforcing electric field. 

It should be pointed out that this effect could be used for the 
measurement of the strength of the quadrupole radiation masked by the 
background in the spectrum of the spectrograph. Then, in favourable 
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conditions (see Rubinowicz 1949), the change in the intensity of the 
enforced dipole radiation due to the interference effect becomes ob- 
servable. 

To detect the interference effect all the experiments should be 
made in the absorption spectrum to avoid the perturbing influence 
of molecular fields. 


The interference effect 


We start from Dirac’s radiation theory. Assuming an atom to 
be in a radiation free space, we have the following expression (see, 
e.g, Bethe 1933, page 429, formulae (38,1), (38, 2)) for the proba- 


> 
bility of a single quantum emission with the wave vector k in a di- 
rection contained in a solid angle dQ and a polarization given by the 
y > 
unit vector e 


> ehy |> >= $2 


-> 
wlk, e) d= S an 6 ‘pe dQ, (1) 
where 
phe te | ikr 
Pyy=<N' |e* 7|ND (2) 


(the labels N and N’ refer to the initial and final states of the atom 
respectively). Since we are calculating the intensities of the Zeeman 
multiplets, we can assume an equal population of the quantum states 
considered and neglect the variability of ». For the intensity I we 
get thus 
>> > 
I(k, e)=const|le Pri. (3) 
(naa _ >> 
We expand now Pyy in a power series in k r and break off with 
the first power expression. A known transformation (see, e. g., Blo- 
chincew 1949, pp. 368, 369) of the two obtained terms performed 
with omission of the magnetic dipole term gives 


>> >> i> 
I(k, e) = const |A Dw +5 Omn) (2, (4) 
where 


> > 
Dyn =<N'|r|N> (5) 


> >>> 
Qnn=<N'\(kr)r|N>. en) 
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i Now we introduce polar coordinates with the origin placed in- 
side the radiating atom (Fig. 1). The angles y, determine the direc- 
tion of observation and the angle a the direction of polarization. The 


Fig. 1 


> 

æ+iy, z coordinates of the vector e can be expressed by the angles 
p, 8, a in the following manner 

€x+iy= (— cos 8 cos a F i sin ajeti”, 7 

ez = Sin 3 COS a. (7) 

. = Ze . 

Expressing Dwy, Qwn by the x+iy, z coordinates, we conclude 
that the only non-vanishing elements are: 


"Dam o for the transition M’=M+1 
D: , Quet- å» č » 5 M'=M (8) 
, Qatiy)? » » j 29° MY’ = M+ 2 


(M= magnetic number). 


This selection rules for Dyy are valid in the case considered, i. e., 
for the electric field parallel to the z-axis. 
Substituting (7) in (4), we obtain with the help of (8) 


I. M'=M+2 


a k? ; 4 = 
I(k, e) =6 5 (Oes ((sin? sin? ats sin? 2 # cos? a} (9) 
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>> i 
I(k, e) = BA{(Dx+iy)? (cos? 0 cos? a+ sin? a) (9) 
+43?(Q(-+ig)2)2 (C08? 28 cos? a+ cos? Ò sin? a)+4kQ(+igyeDxtiy SIN 2a sin? ð} 


Ill. M’=M 


Se Fe 3 k2 l ; 
Ik e) pio sin? ð cos? a + 16 (Qt/(x+ty) (x—iy)—241 )2 gin? 28 sin? at, 


where f is a constant. Thus we see that an interference term occurs, 
but only in the case of AM=-+1 transitions. This term does not change 
the sum of the intensities of both AM =+1 components. The inter- 
ference effect vanishes for observations performed without a nicol, as 
the two quanta with the polarization directions a and —a occuring 
in this case give together an intensity, equal to that without the inter- 
ference term. Henceforth we shall be concerned only with the AM=+1 
transitions. 

The ratio of the absolute value of the interference term to the 
sum of the two remaining terms reaches its maximum in the position 


T D2 


ðr=7; Cos E EO (10, I) 
and if k?/4Q? >D? also in the position 
y D? +4 kQ? 1 k?Q? 
Tamera pagent e9 4 = p o 
For these optimum angles and for M’=M-+1 (9) takes the form 
1 k2 D2 2 
for (ðn an- T=EB p ppU sen(QD) sgn a), (11,1) 
» (mar)... T=4 b D?(1+ sgn (QD) sgn a). (11, IT) 


The simplest way to observe the interference effect is to carry 
out the observations transversaly (ðr). The nicol must be placed first 
at the angle ar and then at the angle —az. The difference of the in- 
tensities in both positions, caused by the interference only, is as fol- 
lows: in the first position one of the two AM=+1 components is 
completely extinguished and the other one twice as strong. In the 
second position both components exchange their roles. One can either 
compare the intensities of the two corresponding components AM=+1 © 
or observe the change in the intensity of one component during the 
change of the nicol position from a; to —a;. The intensity of all the 
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remaining components (A M=+2,0) is the same in both positions. 
It should be noticed that the same effect as by alteration of the nicol 
position (change of sgn a) can also be produced by altering the direc- 
tion of the electric field (change of sgn(QD)). 


Final evaluation of the effect 


We base this evaluation on the following assumptions: 

(1) The electric and magnetic fields are parallel (and parallel to 
the z-axis). . 

(2) The Zeeman splitting is small compared with the term splitting. 

(3) The electric field is weak compared with the magnetic field, 

e., the Stark shift is small compared with the Zeeman shift. This 

means again that we regard the electric field merely as a factor en- 
forcing the dipole radiation and neither splitting nor displacing ap- 
preciably the energy levels. 

4. Russel-Saunders coupling. 

Obviously we can neglect in most cases the perturbation of the 
final state. Then we have (cf. Condon and Shortley 1935, p. 412, for- 


mula (5)) 
> 
<R'L' J’ M'|r|RLJ MY’ 


CR'L'I' M1 PRODI" MCRD Me RET MY (12) 
-r5 


A RL) —H( RD") 
where < >’, < > are the perturbed and unperturbed matrix elements 
respectively. L,J, M, R are respectively the quantum numbers of or- 
bital and total momentum, the magnetic number and all the rest of 
the quantum numbers; F is the intensity of the electric field, s the 
absolute value of the electron charge, H(RL) the energy of the R, L 
state (we neglect. the dependence of H on JM). If we restrict our- 
selves to one electron transitions (where AL= +1 for spontaneous dipole 
transitions) we have non-vanishing terms in (12) only for 
DL+1 for 1’=L 
D=! Dale y DSL- (13) 
E a I ESA 
The order of magnitude of the sum (12) is determined by the term 
with the smallest denominator which is equal to 
E(RL)—G(LZ+1) for L’=L+2, i 
AE=? P(RL)—6(L—1) „ L’'=L—-2, (14) 
E(RL)—6(L+1) „ L'=L, 
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where G is the energy of the nearest state to the RD state with its 
orbital momentum equal in the three cases (14) to D+1, D—1 and 
L-+1 or L—1 respectively. Using hydrogen functions one can prove 
(Bethe 1933, p. 441), that the numerator of the term with the de- 
nominator AE is much greater then the numerators of all the other 
terms occuring in the sum (12). If we assume a sufficient likeness to 
the hydrogen atom, we have the following approximate equation 


-=> 
¢R'L'J' M|r|RLI MY’ 


Po RT J’ M'|r|R'L" I" MD CR'L' I" M\z|RLIM> 


E (15) 


eF IT! T’! ma j 
= <R'LI'M'\r -2|RLIM). 


Using the known formulae for the quadripole matrix elements 
(see, e. g., Rubinowicz 1930) and substituting the expressions (15) in 
the equations (9), we obtain 


1. J —=J—2 
L M =M 


EN I M —any{(< ae 


E 
II. M’=M+1 
ag AA l ef i 
I(k, e)=7 (J? — M?)(J- M —1)(J FM — (5 gy) ator Si 


J sin? 8 cos? a+ = 12 sin? 26 cos? a} By 


keF 
= k?(cos? 20 cos? a+ cos? 3 sin? aye = 2 AE sin 2a sin? al py 


Ill. M’=M+2 (16) 

>> 
1(k,e) = (JM) (JF M—1) (JF M—2)(J=M— 3) 5; (sin? sint a 

1 
AF 7 sin? 23 cos? a) By 
De ead 

I. M'=M 

>> 
I(k, e) = M(J2— -amh T) sinea cost a+ r të sin? 20 cos? a} Bs 
II. eS 


>> 
I(k, ej qT L2M +1) (JF) JM IF a cece 


lis 29 2 2 j 
Ea (cos? 20 cos? a+ cos? e sin? a a) +5 SE TF AAN 2a sin? os, 
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Ill. M’=M+2 
1(k, = (J+ M)(J = M—1)(J-M—2)(J+M+1)— ay sin? sin? a 


Loe 
+ a sin? 28 cos? a) Bo 


Sh Of =l 
I. M’=M 
>> eF 
1(k, e)= (3) (7+ 194 Tupad ea] sin? 9 cos? a 
4+ (3M?—J(J+1))? af sint 20 cos? a} By 
Il. M’=M+1 (16) 
erie | jek 
I(k, e)= 7 CM1)? JFJ EE) (cos? ? cos? a+ sin? a) 
+ = (cos? 20 cos? a+ cos? ® sin? a) + : a sin 2a sin? ah 


Hl WKE 
>> k? 
I(k, e)= (JF M)(JFM—1)(J +M +1) (J+M + 2) (sink? sin? a 
w ; sin? 2% cos? a} Bs, 


where f,,f.,8, are constants and A,,A_,,A — numbers independent 
of M1. 


Estimation of the effect for the first line of the 
? Ds, 7, —>?S:, series in the potasium spectrum 


= The Zeeman effect of the line 2Dy, > 28y, (A= 4642,27-+ 0,02 A), 
2 Ds, > *S1,(A= 4641,77 + 0,02 A) was investigated by Segré and Bak- 
ker (1931), who showed that the line is due to a quadripole transition. 
The intensity of the applied magnetic field was 7,5.10°—17,8.10° Gauss. 
~The terms which are cf interest for us are shown in Fig. 2. According 
to our former considerations we take into account the perturbation 


of the 3D term by the 5P term only. The formula (16) with t=, 


=3 and with f,/f,—6,2.10— (from the experimental results of Segrè 


and Bakker) suggests the application of an external electric field of 
the order of 10 kV/cm as suitable for an observation of the inter- 
ference effect. Fig. 3 shows the Zeeman splitting for this field intensity. 
Our evaluation has been only approximative and we are not able there- 
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fore to fix the exact conditions in which the particular components 
are totally extinguished. We can say, however, that for F of about 


OP Pe : 
=5/9 1347,52 
3d 2D41, 94, 
J=3/ 13470,26 — q 
tp? Py, 1, : 
4s? Sy es he g 


Fig. 2. Term values of KI in cm—1 (from Bacher and Goudsmit 1932 and Grotrian 1928) ; 


> 2 2 
22D, - ‘Sy, Ds, a S ys 
Ep ee 
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. _~ 
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og r a) ) p T A s = er 
a a ish chal ashes WR eR ae b Phi bo Sasi ete. ud Za 
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THE TEMPERATURE DEPENDENCE OF THE WORKING 
OF GEIGER-MULLER COUNTERS 


By Olena STANISZ, Physical Laboratory of the Mining Academy, 
Krakow. 


(received June 26, 1951) 


The measurements were made on selfquenching argon-alcohol Geiger-Miller 
counters in the range of temperatures between —22° and +60°C. The results show 
that the characteristic curves remain normal up to the temperature at which the 
quenching vapour begins to condense, and all changes in these curves can be ex- 
plained by the assumption that the only phenomenon which causes them is the 
condensation of the quenching vapour. 


Introduction 


It would seem that a change in the working conditions of a self- 
quenching Geiger-Miller counter caused by a decrease of temperature 
may occur only at the condensation temperature of the_ quenching 


vapour. Of course the operation of non-selfquenching counters with | 


pure filling gas does not depend on the temperature. However, Korff, 
Spatz and Hilberry (1942) using a counter with the familiar argon- 
alcohol mixture, observed differences in the characteristic curves at 0°C 
as compared with the characteristic curves at +26°C, and a total 
disappearance of the plateau at —22°C. The characteristic curve at 0°C 
is shifted towards lower operating potentials, the counting rate is in- 
creased, and the plateau width shortened. As no data were given on 
the alcohol pressure in the counter and measurements were not made 
in close intervals of temperature, we cannot judge whether the tem- 
perature at which the condensation takes place is really critical for 
the operation of the counter. 

Recently Parkash (1949) and Parkash & Kapur (1950) have been 
investigated an argon-alcohol counter with an alcohol vapour pressure 
of 15 mm Hg in the range between +9°C and +60°C. The authors 
obtained changes of the characteristic curves, i. e., shortening of the 


plateau and high increase in the slope but no shifting of the threshold, 


ied 
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even at +16°C, whereas at +9°C they observed a total disappearance 
of the plateau. They try to explain their results by partial conden- 
sation of the quenching vapour. This interpretation cannot be accepted, 
however, as alcohol vapour at the pressure of 15 mm Hg condenses 
at +2,5°C and in the whole investigated range of temperatures the 
vapour was unsaturated. The fact that the authors do not get any 
shifting of the threshold voltage may be regarded as further evidence 
that the changes in the characteristic curves are not connected, as 
the authors suggest, with condensation of the alcohol vapour. For 
higher temperatures (+58°C) the authors also obtained a shortening 
of the characteristic curves. 

As we see, all the above results on the working of G. M.-counters 
at lower temperatures do not really solve the problem. In addition, 
they involve data which do not agree with the simple assumption 
expressed at the beginning of this paper. 

The present investigation was undertaken in order to find the 
lowest temperature at which the operation of the widely used argon- 


„alcohol counter remains normal, and also to verify the above-men- 


tioned supposition that so long as the quenching vapour is unsaturated 
the characteristic curves of the counter remain constant, 
the changes in them appearing but at the temperature in 
which the condensation takes place. Further, the work of 
the counter in temperatures lower than that of condensa- 
tion was to be investigated, as well as the gradual changes 
in the work of the counter with decreasing temperature. 


Apparatus 


The counters were prepared from brass tubing with 
a wall-thickness of 1 mm, the inner diameter being 2,3 cm. 
Their effective length was 11,5 cm. This is, with a few 
modifications, the type of counter commonly used in the 
Physical Laboratory of the Academy of Mining’. The de- 
tails of construction are shown'in Fig. 1. All the joins, glass | 
jmetal and metal/metal, were made by means of ,,Araldit”,- 
without tin. It should be stressed, that the length 15mm Fig. 1. 
of the blass tube which holds the central wire is important Details of 


for the work of the counters at lower temperatures. ie 
. i F ; ion o: 
The counters were filled with the usual mixture of ma EAS 


G a = 


argon and ethyl alcohol at 90 mm and 10 mm of mercury 
respectively. The temperature at which the alcohol vapour condenses, 
is then —2,5°C. 

1M. Miesowicz et al., unpublished. 
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Each of the counters ended in a glass tube fixed into the collar 
of the counter with rubber. The joins were strengthened with isolating 
tape and then coated with shellac, so that when the counters were 
immersed in cooling mixtures the anode was effectively isolated. The 
counters were connected in groups of three and in the midst of each 
group a weak radium preparation was placed. 

The counters were cooled by means of eutectic mixtures. This 
mode of cooling ensures constancy of temperature for several hours 
over the whole space occupied by the counter. The counters were 
completely surrounded by the eutectic mixture which was contained 
in a vessel isolated from outside by a thick layer of sawdust. 

The counter was supplied by a stabilized H. T. supply. The 
pulses were registered in a scale of 16 and 32 on a Baldinger scaler. 

The maximum was 100 pulses per second, so 
BN that corrections of the number of pulses regi- 
stered for the dead time of the counter could 
be disregarded. 


Results 


The measurements were carried out with 
7 counters in three series; in the 3-rd serie, 
besides the new counter two previously in- 
vestigated counters were used as controls. The 
measurements were taken within a tempera- 
ture range from —21°C to +60°C> 
Each characteristic curve was taken twice 
as a check of the stability of the counting 
rate curve. In general two such curves were iden- 
600. so 700 750 #7 tical. For lower temperatures, however, where 
i condensation takes place, it happens that two 
S aa vraat such measurements result in different characte- 
RE A E ARA ristic curves. Then the second curve lies as 
operating voltage, (b) af- ® rule lower than the first. After the counter 
ter some minutes of wor- has been working for some time at a fairly 
King of the counter high voltage above the threshold potential,the 
(c) stable characteristic, counting rate becomes stable. Fig. 2 illustrates 
mes the above phenomenon for a curve obtained 
at —21°C. The dotted curves a,b are obtained immediately after the 
switching on of the operating ove The continuous curve ¢ is the ~ 
stable one obtained after the counter has been working for some time 
of the order of some minutes. 


Fig. 2. Characteristic cur- 
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The results obtained in these experiments were similar in all 
the series of measurements. The results for one counter, which are 
thus typical for all the counters, are shown in Fig. 3. The numbers 
beside the curves give the succesion-in-time of the measurements and 
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600 650 700 750 800 850 900 950 14000 4050 #00 1150 


Fig. 3. Typical diagramm of temperature dependence of the Geiger-Miiller counters. 


indicate. that each measurement at a low temperature was checked 
by a measurement at room temperature. . 

The main results shown by these curves are as follows: 

(1) The counter operates normally after remaining for some time 
in low temperatures, even while working; its counting rate curve at 
room temperature returns to normal. 

(2) All the changes in the course of the curve appear only at the 
temperature at which the quenching vapour begins to condense, i. e., 
in the present case at —2,5°C. In the range of temperatures from 0° 
to +60°C complete agreement in the curves is obtained. 

(3) The following changes in the characteristic curves were ob- 
served: 

(a) A shifting-of the threshold towards lower potentials. 

(b) A shortening of the plateau of the curves. 

(c) An increase of the slope (in the case of some measurements). 

(d) A raising of the plateau. 

(4) The counter can work in a lower temperature than that of 
condensation of the quenching vapour. A fairly long plateau in the 
characteristic curves may be obtained even below this temperature. 
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The fact that the changes in the characteristic curves begin at 
the temperature in which the quenching vapour begins to condense, 
suggests the following interpretation. Cooling of the counter not below 
the temperature of condensation of alcohol vapour does not cause any 
visible effect, as the change of pressure under the influence of de- 
creasing temperature has only a negligible effect on the free path of 
the electron bringing about the avalanche so long as the number of 
molecules remains constant. This is also why counters filled with pure 
gas show no temperature effect. 

The shifting of the treshold towards lower potentials can be ex- 
plained by the decrease of the amount of alcohol vapour caused by 
condensation. When the number of molecules begins to fall, the free 
path of the electron bringing about the avalanche increases and the 
threshold voltage decreases. The same effect is noticed when the per- 
centage of alcohol vapour in the counter is reduced. 

The decrease in the number of molecules of alcohol vapour ac- 
counts also for the increase of the slope and the shortening of the counting 
rate curve. When the number of alcoho] molecules decreases the num- 
ber of argon ions reaching the cathode increases. It is known that 
this causes a rise in the rate of the multiple and spurious pulses and 
as a result the characteristic curve becomes shorter and steeper. The 
shortening of the plateau and the increase of the slope may be also due 
to the appearance of conducting domains of condensed vapour bet- 
ween the' electrodes which cause additional discharges to be possible. 
The decrease in the number of alcohol molecules caused “by conden- 
sation is a reproductive process; on the other hand, the process of 
discharge through liquid is much harder to seize on account of the 
haphazard way in which the conducting domains arise. 

It seems that the appearance of the raised characteristic curves 
at temperatures at which alcohol vapour is saturated is connected 
with the arising of additional pulses caused by the presence of liquid 
in the counter. Observation on the oscilloscope has shown that below 
the temperature of —2,5°C their number increases with decreasing 
temperature. These pulses are small and only some of them are re- 
gistered by the apparatus. The upward shifting of the plateau ob- 
served for lower temperatures would then give the number of ad- 
ditional pulses registered (Fig. 3, curves 14 and 8). z 
; The results of the above measurements show that all changes 
A the characteristic curves can be explained by the assumption that 

e only phenomenon which causes them is the condensation of the 
quenching vapour. In contrast to the results obtained by Parkash 
the changes in the characteristic curves always appear at the egm 


Be aia B 
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perature of condensation and are always accompanied by a shifting 
of the threshold which is a sure proof that the process of condensation 
has begun. There seems to be no contradictions in such a picture. 
I want to express my sincere gratitude to Professor M. Miesowicz 
for suggesting the subject and for helpful discussion of my results. 
My thanks are also due to the Sekcja Naukowa Komisji Popierania 
Twórczości Naukowej i Artystycznej przy Prezydium Rady Ministrów 
R. P. for financial assistance. 
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INFLUENCE OF THE ELECTROLYTIC DISSOCIATION 
AND THE HYDRATATION OF H,SO,— MOLECULES ON 
THE SCATTERING OF SLOW NEUTRONS 


By J. A. JANIK, 2nd Physical Laboratory, Jagellonian University, 
Krakow 


(received July 6, 1951) 


Coefficients of absorption of slow neutrons in H,SO,— water — solutions 
were measured. At low concentrations of acid the effect of electrolytic dissociaton 
was observed. The results are in agreement with the theory if one assumes com- 


plete or almost complete dissociaton in H+. and SOF ions up to the concentration 
of 40 per cent (by wt.). At higher concentrations the increase of the absorption 
coefficient due to the hydratation of H,SO,— molecules was obserwed. 


1. Introduction 


It is well known that the scattering of neutrons by protons de- 
pends on the chemical binding of the protons in the molecule. Acco- 
rding to Fermi’s theory (1936) the cross-section for scattering of zero- 
energy neutrons by protons elastically bound with an infinitely great 
mass is four times as great as for free protons, according to the formula 


Dy GO? Moo 


where oy is the cross-section for the scattering of neutrons by free 
protons, dœ — cross-section for the scattering of neutrons by protons 
bound with an infinitely great mass, u — reduced mass of the neutron 
plus molecule system. 

This theory is valid for zero-energy neutrons, i. e., for neutrons 
of much lower energy than the vibrational and rotational 
the molecule. 

A proton may be considered as free when the energy of the in- 
cident neutron is large in comparison to the rotational and vibrational 
quanta. The cross-section for the scattering of such neutrons by pro- 
tons is (Cohen et al. 1939; Rainwater et al. 1948, Jones 1948) 


on=20-10-%4 cm2. 
| R d 


quanta of 
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In the intermediate case for which the energy of the neutron is 
great in comparison with the rotational but low in comparison with 
the vibrational quantum (Hyot<E,<Eyn; En >Erot>kT"), the cross- 
section for slow neutrons can be theoretically obtained from to the 
theory of Sachs and Teller (1941): 


~ 214 
T= Col fy mu S parE nib, | , 


where Mı, Mo, are the characteristic values of the tensor defined by 
pti=14+mM=14+%n 


M is the mass tensor of a molecule, m — the mass of the neutron, 


o — the length of the vector m? Pi, P; — coordinates of the momen- 
tum of the molecule, f; — its direction cosines, a, — direction cosines 
of the momentum of the incident neutron. 

It may be seen from this formula that the neutron cross-section 
depends on the structure of the scattering molecule, on the direction 
of the neutron beam, and on the neutfron energy. In the case when 
we are interested in the total cross-section for the scattering of ther- 
mal neutrons, we are obliged to take an average for all directions 
a, and fi, paying due regard Maxwell’s velocity distribution of these 
neutrons. 

At last, when we can neglect the influence of the rotation but 
when the neutrons are not of zero-energy, the scattering cross-section 
is given according to Bethe (1937) by the formula 


CH= Too 


L EAN 1 
ET (Eate) =], 
where : 
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E, is the energy of the neutron and @,,@,. are the vibration frequencies 
of the proton in the molecule. 

All these theoretical considerations were confirmed experimen- 
tally recent years. In these experiments chemical bonds, intermole- 


cular forces, and similar molecular problems were investigated by the 


help of the scattering of slow neutrons by bound protons. — 
It seems to be interesting to investigate in this way the con- 
ditions which exist in the solutions of sulphuric acid of various con- 


centrations. The idea of this problem is the following: One measures 


the absorption coefficient due to the scattering of thermal neutrons 
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in water and in pure H,SO,. If in the solutions of sulphuric acid, only 
H,O and H,SO,-molecules were present, it would be possible to cal- 
culate the absorption coefficient for every concentration and to com- 
pare the obtained results with the experimental data. But we may 
suppose that in the solutions of sulphuric acid not only H,O and H,SO,- 
molecules are present but also HF, SO; and HSO;z-ions as well as 
molecules of the hydrates of sulphuric acid of the type H,SO,.nH,0. 

One could expect that in H,SO,-solutions, especially at lower 
acid concentrations, the effect of the presence of the Ht-ions could 
be experimentally verified, as the scattering cross-section for a free 
proton is smaller than for a bound one. A positive or a negative result 
of the measurements had to decide whether the protons of the disso- 
ciated electrolyte can be considered as free or not. 

It was clear from the outset that only strong acids could give 
an appreciable effect as these acids are completely or almost com- 
pletely dissociated already at great concentrations. But it was not 
possible here to use any strong acid with a complicated molecule, 
since the molecular mass could not be very great as this would pro- 
duce too small a difference in the scattering cross-sections for dis- 
sociated and undissociated molecules. 

All these considerations and the facility of obtaining sulphuric 
acid in graet concentrations decided for its choice as the neutron scat- 
terer in these experiments. 


2. H.SO.-s olutions 


The structure of the solutions of sulphuric acid is not yet known 
quite well. Investigations of many authors give a rather qualitative 
than quantitative picture of this structure. One of the most important — 
is the work of Cartwright (1936) on the absorbing and reflecting power 
of H,SO,-solutions in the far infrared. Another important work are 
the investigations of Woodward and Horner (1934) of the Raman 
spectrum of sulphuric acid solutions at various concentrations. 

Both these works give us the following picture of H,SO,-solutions: 
the H,SO,-molecules added to water are mainly dissociated in two 
Ht-ions and a SO; ~-ion untill the concentration about 12n (43 per 
cent by weight) is reached. With increasing concentration the ions Ht 
and HSO; begin to predominate in accordance with the above men- 
tioned Raman spectrum investigations. Both these works do not give 


us any quantitative suggestion concerning the degree of dissociation 
for the first and the second kind of ions. i 
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Another important property of H,SO,-solutions is the building 
of hydrates. The dependence of the freezing point of H,SO,-solutions 
on the acid concentration show several maxima which are interpreted 
as the influence of the hydrates: H,SO,.4H,O, H,SO,.2H,O, and 
H,SO,.H,O. 


3. Apparatus 


In this work the slow neutrons were counted by means of a boron 
lined ionization chamber. The whole apparatus was described in a pre- 
vious work on the absorption of slow neutrons (Janik 1951). 

Fig. 1 representing the electrometric stage of the amplifier with 
the ionization chamber shows the general geometry of the means- 


scatterev 


3V 
Fig. 1. The electrometric stage of the amplifier with the ionization chamber. Dis- 
tances in cm. 


rements. A source of photoneutrons slowed down by 6 cm of paraffine 
was used. The necessary y-ray source consisted of 200 mg of radium 
enclosed in platinum tubes and needles. The scatterers were placed 
in rectangular glass vessels which were situated between the chamber 
and the neutron source in a distance of 11 cm from the ionization 
chamber. 


4. Results 


The absorption coefficients for scattering of slow neutrons by 
—H,SO,-solutions of 96, 70, 60, 50, 40, and 32 per cent (by wt.) con- 
centration and by destilled water were measured. In order to avoid 
electromagnetic and mechanical disturbances the measurements were 
taken only during the night. 


The results are represented in Table 1 and Figs 2 and 3. 
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Per cent of 
H,S0, 
in water 
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Table 1. 


(cm!) 
measured. 


2,726+ 0,145 
2,326 + 0,060 
2,245 + 0,062 
2,289 + 0,076 
2,489 + 0,064 
2,346 + 0,094 
1,411+ 0,067 


u (cm—!) 
caleulated for 
undissociated 
H,S0,-solutions 


2,726 + 0,145 
2,566 + 0,139 
2,483-+ 0,135 
2,385 + 0,133 
2,256 + 0,126 
2,096 + 0,120 
1,411 + 0,067 
1,261 + 0,080 


The measured absorption coefficients given in the second colu mn 
are also plotted in Figs 2 and 3, curves A. The absorption coefficients 
given in the third column were calculated from the measurements 
for H,O and a 96 per cent H,SO,-solution which may be considered 
as practically undissociated (Figs 2 and 3, curves C). w 

If we denote by u and üa the T coefficients (in cm- a 
of the constituents of a mixture, by o, and o, their densities, by ` 
and p, their aboundances in the mixture, and by ọ the density of the 
mixture, the absorption coefficient u of this mixture is given py the 
well-known formula > 
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The theoretically calculated cross-section for a water molecule on 


the assumption that the vibrational frequencies of a proton in the H,O- 


molecule are w= 3600 cm— and w,=1600 em— is ono =126,7.10—*4 cm. 


This theoretical value is much greater than that obtained experimen- 


25 
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Fig. 2. Dependence of the linear absorption coefficient on the concentration of H,SO,; 
(per cents by weight) A. Experimental results. B. Theoretical comparison C. Eva- 
luated from the mixture rule. 


tally. Rossel suggested that the reason for this discrepancy is the 
neglecting of the rotation effect of water molecules. Taking this effect 
into account, we obtain, MEO OY to the formula of Sachs and Tel- 
ler (1941), 


*~ 


0n,0= 78,5. 10—*4 cm?, 


which stands in a satisfactory agreement with the experimental results. 
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It may be seen from Table 1 and Figs 2 and 3 that the figures 
given in column 2 of the table, calculated on the assumption that 
only H,O and H,SO,-molecules are present in the H,SO,-solutions do 
not coincide with the measured absorption coefficients. For the con- 
centration below 40 per cent of H,SO, the measured coefficients are 
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ions can be considered as free protons with a Scattering cross-section 
for slow neutrons of 20.10-*4 cm? only. Therefore, the scattering cross- 
section for the sum of the H® and SO -ions ii be (20+ 20+1,5+ 
+4.3,5).10—-*4em?= 55,5.10—*4 em? (os=1,5. 10—24 em?, co= 3,5. 10—44 cm? 
Goldsmith et al., 1947), which corresponds to. the absorption coeffi- 
cient u= 0,626 cm—. It is 2,02 times smaller than for the same eon- 
stituents bound in the H,SO,-molecule. The calculated coefficients on the 

MHSO, 
(HF +HT +80) ~ 992 
are plotted in the curves B of Figs 2 and 3 (below 40 per cent). There 
is a satisfactory agreement with the coefficients obtained in the meas- 
urements. 

With increasing concentration of the acid the influence of the 
dissociation becomes smaller, but in the same time the hydrates 
H,SO,.4H,O, H,SO,.2H,O, and H,SO,.H,O begins to appear. The 
associated water molecule has a restricted possibility of rotation which 
causes, according to the theory of Sachs and Teller, an increase of 
the scattering cross-section. If we put o, o= 126,7. 10—*4 em?, (tu, = 
= 4,215 cm—"), we obtain for the H „SO,-solutions of 60 and 70 per cent 
the ao plotted in the curves B of Figs 2 and 3 above 60 per cent. 
These values are much greater than the experimentally obtained ones. 
Therefore it may be supposed that the association of the water mole- 
cule with the H,SO,-molecule does not make the rotation impossible 
but only renders it difficult. This may give some evidence of 
a very loose chemical binding between H,SO, and the associated 
H,O-molecule. 

I. wish to thank Professor H. N iewodniczanski for interesting and 
helpful discussions, as well as Professor W. Kemula, Professor K. Gu- . 
miński, and Mr. K. Pigon for many informations concerning the struc- 
‘ture of the H,SO,-solutions. I wish also to express my gratitude to 
Dr. E. Wyrobek for his courtesy in permitting me to use 200 mg Ra 
belonging to the Instytut Onkologiczny in Cracow. 


assumption that Myo = 2,726 em and u = 0,626 cm—! 
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FIELDS DEFINED BY ELEMENTARY LAWS 


By A. RUBINOWICZ, State Institute of Mathematics, Warsaw 


(received July 29, 1951) 


Supposing that a given scalar field u is isotropic and homogeneous and cor- 
- responds to a linear physical law it is possible to define it by an elementary law (pos- 
sibly empiric) in the form of a spherically symmetric field of a simple source u= f(r) 
Its partial fields are namely given by f(r) and the derivatives of f(r) with regard 
- ġo Cartesian coordinates x,y,z. The fields of the l’-th derivatives may be interpreted 
as fields of multiple sources of order 2. 
By a linear superposition of such fields it is possible to obtain peels 
fields F; m(0, p)Ri(r ) of order 2’ (where V =1, 14 2, 1 + 4...) given by simple formulae + 
and having generally singularities in the source point r= 0. p 
To obtain functions of the form Y;m(0 É) Si (r) which are finite at r= 0, 
the field of m simple source u= f(r) may be dëveloped into a series containing the 
functions Ki (r). 

In the field given by a superposition of all these fields Y,,(6,p) R}(r) 
and Yim (6.9) Si (r ) it is possible to formulate boundary value problems for 
a sphere, if the number of linearly independent functions Ri (r) and S$ (r) with a gi- 
ven J is finite. 

We may AEA put the question under which suppositions it is possible 
to formulate boundary value (and in certain cases also eigenvalue) problems for 
domains of an arbitrary shape for such a field given by an elementary law only. 

The formulae expressing the multipole fields Y; m(0 p) Bi (r) by the fields 
of multiple sources! can be applied also to write down immediately all the multi- 
pole solutions and in certain cases also all the solutions regular at r=0 of any dif- 


n 
ferential equation of the form 3X asAs u=0 (especially the potential or the wave 


equation). The solutions winch arc regular at r=0 are obtainable in any case also 
in another way (§ 5). It is further possible to use the above formulae to replace 
in Green’s formula a simple source by arbitrary multipole sources using simple de- 
rivations only. All formulae apply also to time-dependent fields, for example in 


the case of the wave equation. 
In an appendix the question was discussed how far a given field determines 


its source system. 


§ 1. Definitions ¿and general remarks. The aim of this 
paper is to discuss the characteristic properties of linear scalar fields 
defined in an isotropic and homogeneous space by the field of a sphe- 


1 Cf. (13) and (18). 
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rically symmetric simple source or by an elementary law, possibly given 
empirically. We suppose therefore in the following considerations that 
we have to do with a scalar field u with the following properties: 

(a) There exists a spherically symmetric field u= f(r), possibly 
depending on the time, which has such a singularity at r=0 that it 
can be interpreted as the field of a simple source at that point. In all 
other points of space f(r) is finite. 

(b) The form of the function f(r) is independent of the position 
of the source point in space. 

(c) The principle of superposition applies to the given field. 

Under these suppositions the l’-th derivative 


TE OBI) as 
U= Jya gy? dee” V’=a+b+e (1) 


represents also a field which has in general a singularity at r=0 and 
can be interpreted as a field of a multiple source of order 2” generated 
by 2” positive and negative simple sources arranged in a definite man- 
ner. We obtain namely a multiple source of order 2” if we take a multi- 
ple source of the next lower order 2”—! and place opposite it a multiple 
source of the same order obtained by a translation of the first source 
and a change in the sign of its strength. There exist 


g wel ae) A 


such multiple sources of order 2” (cf. Appendix). 

We call the variety of fields we get by a linear superposition of 
the fields (1) simply the fieid belonging to the field of a given simple 
source or a fundamental law u= f(r). But, as we shall see later (§ 5), 
in some cases we shall have to supplement these fields belonging to 
multiple point sources by fields generated by spatial or surface distri- 
butions of simple or multiple sources which are generally finite at r=0. 

It must be remarked, however, that it is also possible, of course, 
to build up fields (1) with functions which are finite at r=0, but in 
this case higher derivatives (1) may become infinite. 

The notion of a multiple source is associated with an intuitive 
geometrical and physical picture. But to obtain an easy survey of 
the symmetry properties connected with the isotropy of the field (1), it 
is advantageous to introduce a more abstract notion namely that of 
a multipole source of order 24. We define its field as being such a linear 
combination of the fields (1) of multiple sources of order 2” which 
has the form j 


u =F (0,9) Ri(r). (3) 
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Here Yım(0,p) is a spherical surface harmonic. The radial function 
Ri (r) does not depend on m and has for r=0 generally a singularity; 
it has two indices J and 7’ to indicate the order of the multipole field 
and the order of the fields of the multiple sources of which it is build up. 
The number of the multipole fields of order 2! (belonging to a given 1’) 
is 21-1. These 2/+1 fields are all linearly independent. 

All fields of multiple sources of the same order 2” change or do 
not change their sings under an inversion with respect to the source 
point according as l’ is odd or even. Just so a multipole field (3) chan- 
ges its sign or remains unchanged under an inversion according as 1 is 
odd or even. It follows from these facts that both indices J and 1’ in (3) 
are of the same parity. 

From (2) we conclude that Ay—Ay_»2 is equal to 2I’+1, i. e., to 
the number of the multipole fields (3) of order 2”. Ay is therefore equal 
to the sum of the numbers of multipole sources of the orders 2”, 2-2, 
2—4... This fact leads us to assume that by a suitable linear combi- 
nation of the fields (1) of the multiple sources of order 2” we can ex- 
press any multipole field (3) with J=I1’,l’—2,l’—4,... Therefore, the 
order 2! of the multipole fields should not be greater than the order 2” 
of the multiple sources of which it is build up. In (3) we have then 
always 1<l’ and because Z and l’ are of the same parity all Rl(r) be- 
longing to the same spherical harmonic are given by 


Ri *(r), (s =0,1,2,...). (4) 


We have therefore generally an infinite number of radial functions (4) 
belonging to the same spherical harmonic Y,,,(6,¢). 

From the point of view of group theory the dy fields (1) form 
a Ay-dimensional representation of the rotation group which is re- 
ducible into the 1’, l’—2,l’—4,...-dimensional irreducible representations 
of this group given by the multipole fields (3). 

We define the spherical surface harmonics in (3) by 


Yin(9,g)=emePT(8), aa (5) 


where the associated Legendre functions P7'(6) are not given as de- 
fined by Franz Neumann but by C. G. Darwin 


m (l m)! +m qi+m 2 I 
P (6) =- sin 9 Teos 6m gji (008 0—1) | 
! ž i —m)(l—m—1 


(t—m)(1—m —1) (l—m —2) (I— m— 3) ET AE 
i TETE Ee 7 ) 
A 11* 
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Speaking of an elementary law in physics we usually think of 
two infinitesimal parts of one or two different physical systems exerting 
an influence on each other, e. g., an infinitesimal current element and 
a magnetic pole or two equal or different elementary particles. But 
in the following considerations we shall use this word also in 
a broader sense, e. g., to denote the fundamental solution of a diffe- 
rential equation as the potential equation describing the velocity 
potential of a fluid flow. 


§ 2. Multipole fields expressed by fields of multiple 
sources. Our problem is to find out for every multipole field (3) of 
order 2! a linear superposition of fields (1) of multiple sources of or- 
der 2” having the directional symmetry of the field (3). The right tool 
to do it is the group theory. 

First we try to express the multipole fields (3) of order 2? by 
the fields (1) of multiple sources of the same order, i. e., for J=l’. We 
denote for that purpose by a,b,c the components of a zero vector, 
so that 


+ B+ &=0. (7) 
Taking into account that (ax pos tos pee (ax + by + d 
on oy \ de) or Yat CU 


ə ə ə 
that C (ax + by + cz) = a + b+ c=0, we get the 


relation : 
ə ə ə \! IE GNE 
(azz tta teg) = awt dy + ee) (=F) fr). (8) 
It follows from (7) that the homogeneous polynomial in æ,y,z of degree 1 


(ax + by + ez) (9) 


is a solution of the potential equation and is therefore equal to a pro- 
duct of r' and a linear combination of spherical surface harmonics 
of order l. : 

We can fulfil equation (7) identically by expressing a,b,c by 
means of the spinor components u,v: i ki 


u2— y2 EN y2 


a= y 55 c= — we. 


In this way (9) becomes a homogeneous polynomial in u,v of degree 21. 
After introducing here instead of w,y,z spherical coordinates r,0,p by 


a+ ty=ret? sin 6, z=r cos 6 (10) 
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the coefficients of Weyl’s monomials w= y+ in (9) become spherical 
harmonics. We get namely 
l 
(ax by + cz)!= (—1)! r! N ym ypttm l! 
T = a (L—m)! lm)! ¥i,m(9;9). (11) 


m=—l 


Here we can put according to (5), (6) and (10) 


(21)! aeaaea Atel Mm —1) 
WY 1m 9,9) = ory (x+ ty) fa TA 2(21—1) Pai i 
TAE ON m—1)(l—m—2)(l—m—3) ae a) 
a ledel 3) n =: 
By the same procedure the operator 
ð ð ð 
(attat 35] 
goes over into an operator we get by inserting in (11) expression (12) 


in which we have substituted for x,y,z and r? the operators x = 


and 4. From the equality of the coefficients of the monomials um ym 
on both sides of (8) we get finally the formula we intended to obtain 


1 d\! 
TAA ie DAS 
Yim, p)" ES) 
ao O FOS _ (=m) (t—m—1) , dom 
ett! zeta ogi a (2T) omaa 


((—m)(l1—m—1)(1—m—?2)(l—m—3) ,, 3m 
oo OA OA =a, 1 ee Air). 


(13) 


+ 


The right-hand side of this equation consists of a sum of fields 
of multiple sources of order 2". Its left-hand side has the form of a multi- 
pole field (3) of the same order with the radial function 


Ri(r)=r'(= z) f(r). (14) 


(13) gives us therefore an expansion of a multipole field of order 2! 
into fields of multiple sources of the same order. 

Formula (13) is applicable for positive m’s only. For negative m’s 
we must use in (5) the equation Pr (0)=(—1)"P7'(8) (which applies 
to the associated Legendre functions in Darwins seeder so that 
in (13) we have to aga the factor (—1)™ and ae m by |m] and t 
by —i. 
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The relation (13) was obtained in a different way by Hobson (1931, 
p. 175). Nevertheless I have reproduced here the way in which I have 
found it because it seems that it can be more easily generalized to 
non-sealar fields. Hobson writes equation (13) in the simple form 


Paleo ee ae] 
i, 952)2"( Fe] 0) = Vil Se Far Ze | HO 


where Uj(#,y,z) is any solid spherical harmonic of order l. But he 
does not make any application of this relation. 

To get multipole fields (3) with 1’=I-+ 2s let us now suppose 
that we have expressed the multipole field (3) of order 2! by a sum 
of fields of multiple sources (1) of order 2” and let us apply the ope- 
rator A to this relation. The application of A to a multipole field (3) 
of order 2! yields a multipole field of the same order 2! 


=. > ak ae dae NERE 
AY uml 9,0) 81 (0) = ilb (5 (2-5 RE (r)) HS? o). as) 


But by applying A to a sum of fields of multiple sources of order 2” 
we get fields of the same kind of order 2"+?. It is therefore possible to 
define 


AY im( 9,9); (7) = Yim(9,¢) Rit (r), 


where according to (15) 


tye Ving! d oy, WLI) ap fet 0 Le ay, 
RP= aghe g EO) Rn) =r ARE) T 
1 
od NEE cd bo SS 
Sas ern ( R! (r)). 
The operator A applied to a multipole field (3) does not change its 


order 2! but raises the index J’ in the radial function R}(r) by two, 
so that 


A° Yim 0,9) Ri (r) =Yim( 9,0) Ry (r). ALD 
Applying therefore A to equation (13), we obtain 


Yim) = Yam Oye) (20+ 3)r (FE) poh B gy 


r ar 
=e ix)" gm  (l—m)(l—m—1) | 3m- 
201 \ðu " dy} Ozma 9 (21—11) A ggat tAn) 


as the operators on the right-hand side of (13) commute with A. 
We Bos therefore the multipole field of order 2! with the radial 
function Ri (r) and the spherical harmonic Y im(9,¢~), if we use in (13) 


Fito 
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instead of f(r) the radial function Af(r). This procedure may be con- 
tinued. To obtain the eres fields of order 2! palin the spherical 
harmonic Yia 0500): Yor (7); Ying tiy W) nik Sp Jie we have to 
apply the operator appearing on the right-hand side of (13) to Af(r) 
A*f(r), A’f(r),... From this fact it follows that also on the left-hand 
side of (13) we can always unite the operator A with f(r). We get therefore 
1d Imai G o GNE 

r £) A* f(r) =r =| (= AR a 5 f(r). (18) 
(18) yields generally an infinite number of multipole fields (3) of the 
same order 2! and with the same spherical harmonic. 

In such a way we have found all the radial functions Rr) 
(s=0,1,2...) belonging to any spherical harmonic Y,m(0,p) of order 2, 
and a given elementary law f(r). Their number is generally infinite, 
as we have suspected, unless f(r) is a solution of one of the differential 
equations A‘’w=0 with a finite index 7 or ‘there exists between the 


RM (1) = 


A*f(r) a linear dependence of the form Sa, As f(r) <0. 
s=0 


If in formulae (18) we lower the index l by 2,4,6..., we get all 
the multipole fields of the orders 2!,2'-2,2!-4,2-6... which are expres- 
sible by the fields of multiple sources of order 2. 

It was the mathematical, especially the group theoretical, aspect 
of the problem which led us from the fields of multiple sources (1) to 
the multipole fields with the radial functions (18). But we can also 
start from a physical point of view. In the case of quantum theoretical 
fields we can look for linear combinations of the fields of the multiple 
sources (1) which are eigenfunctions of the operators of the square 
of the total angular momentum and its z-component. It seems possible, 
for example, to treat in such a manner the quantum theoretical spinor 
fields. 


§ 3. The cases of the potential and the wave equations. 
We apply now relations (13) and (18) to fields satisfying the potential 
or the wave equations Au =0 or Au+ k'u =0. 


In the case of the potential equation we have f(r)=1/r and hence 
by (14) 
1d 1(21)! 1 
Ri(r) En ($ E r = (1) OY yi : (19) 


A(1/r) =0 all terms except the first vanish on the right-hand side 
of (13) and we get therefore in au case of m>0 the expression 


Yiml 9,9) : oa ia 
terre it), (= + di Dzmm y` 
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However, this well-known representation of the potential multipoles 
is not unique. Since 1/r satisfies the potential equation, we can re- 


ð 9 O o Bae 
place here zz JA by -5 Stig] & ~iz| and get so the more gene 


> 


ral formula Pe 
aA  .d\t*/9 Awo Yim(9,¢) i 
(5: dy) \ox dy] dozer r o 


which is applicable also in the case m<0. The non-appearance of any 
constants besides (—1)4+¢ is due to the fact that we are using spherical 
harmonics as defined by C. G. Darwin. 

All multipole fields which can be derived from (13) by continued § 


application of the operator 4 vanish by (18), as f(r)= =i is a Botid 


of the potential equation Au= 0. To a given spherical harmonic Y, sm 9,9) F 
belongs therefore only one multipole field, namely the one with Ryr). 
All other Ri(r) with l'l vanish. 

In the case of the wave equation it is advantageous to introduce 
instead of ©,Y,% and r the variables 


E=kn,. n= ky, C=ke and o 


Using now 4 as short for =z A) = ee = T 5 an the wave a gets the 


form 4u+u=0 and we can use formulae (13) and (18) form erly 
obtained replacing in them #,y,2 and r by ¢,n,¢ and o. For the fiel d 


of a simple source we can take now EA jä 
k . LEE | 


as e o ape 


~ Denot noting as usually — yey 3 ai te aaa n sch OE, 
cat ats reared ala ni 
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On the right-hand side of (13) we have to make use of A(e ®/o) 
=—e—e/o. This means that we have to put —1 for A. Hence e m>0 


a) 


Ven = Yim( O) 


a UE (lL—m) (l—m—1) 3m- 


Sss ce a) bee Ss, 
a a oy) \acen 2 QII) gm 
pliant m1 )(t—m—2)(l—m—3) 3m4 Tene. 
Pet (OI==1) (21-3) acm jeer 


Using (18) we obtain always formula (22) again. Mym and e—le/o 
are namely solutions of the wave equation 4u——u, so that the appli- 
cation of A to (22) causes only a change of sign on both sides of this 
equation. This means that 


Ri (o) =(—1)* Bile). (23) 


By (21) we get in the case of the wave equation the multipole 
fields (22) with a singularity at ọo=0. To obtain regular fields at ọ=0 
we have to use instead of (21) the real part of it, namely f(e)=sinog/oe. 
We get so instead of (22) an expression in which the function ¢,(0) 
is replaced by its real part yo) (cf. (42)) which is regular in o=0. 
On this example we see that it-is possible to obtain by (13) and (18) 
also fields which are regular at ọ=0. 

In the case of the differential equation 4u+ k?u=0 where k is 

an aes or a complex number, we have to use the elementary 
laws f(r) = et"”"/r and f(r) =(e#"—e'*r)/r to obtain by an application of 
(13) fields corresponding to multipole sources at r=0O or being 
regular at that point respectively. 
i But there are also regular fields of the form Y,n(6,¢) Si(r ) be- 
longing to certain differential equations which are not obtainable by 
the method used so far. As an example we quote the solutions Y,,(9,~)r 
in the case of the potential equation. We obtain them in the limit k—0 
from the solutions of the wave equation (cf. equation (42)) which are 
regular at r=0. Now, 


lim #1. 3... (21+ D Vil 0,9) A = Grane (24) 


This solution of the potential equation yields for /=0 the spherically 
symmetric solution f(r)=1= const. If we apply to it the formulae 
of § 2 we obtain vanishing fields only. Nevertheless we may consider 
to a certain extent this spherically symmetric solution as the elemen- 
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tary law belonging to the solutions (24). In §5 we shall show how 
to obtain the solutions (24) starting from the elementary law dir: 

In the cases treated above the reduction of the infinite number 
of multipole fields (3) of order 2! to one multipole field of this order 
is caused by the fact that these fields are solutions of the differential 
equation of the second order Jdu=0 or Au+ ku=0. 


§ 4. Fields belonging to differential equations of higher 
orders. We have always to expect similar reductions of the infinite 
number of multipole fields of a given order 2! to finite numbers when 
the given field is a solution of a homogeneous linear differential equa- 
tion with constant coefficients of the form 


(Xa, A*)u=0. (25) 
s=0 
If we apply namely the operator of this differential equation to its 
solution Yim(0,~) Ri (1), we get by (17) in any case the following re- 
lation 


Z aski (7) =0 


between the n+1 radial functions R}+™(r) (s=0,1,2,...,n) belonging 
to the same spherical harmonic Y,,,(6,~). It is therefore possible to 
express all non vanishing radial functions Rj**(r) by n radial func- 
tions, e. g., by the first n radial functions with s=0,1,.5,n—1. 

But n is only the maximum number of the multipole fields (3) 
belonging to a given simple source f(r). Generally, however, the dif- 
ferential equation (25) will have several linearly independent spherically 
symmetric solutions which can be used as fields of simple sources f(r). 
To determine the number of all the linearly independent multipole 
fields with the same spherical harmonic we must therefore look for 
all spherically symmetric solutions of (25). 

For this purpose we write equation (25) in the form 


an [| (A+k)u=9, 
j= 


where the generally complex numbers kj are given by the roots sj=ki 
of the algebraic equation 


2 a{—#) =0 (26) 
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with the same coefficients a, as in the differential equation (25). If we 
denote by u; any solution of the differential equation 


A Uj + ki Uj= 0; 
then 


n 
u= X uj 
j 


is certainly a solution of (25). A spherically symmetric solution of (25) 
is therefore given by 


o= S (oF +4) (27) 


if none of the k; vanishes. This expression is the general spherically 
Symmetric solution of (25) as the number of the arbitrary constants 
in (27) is equal to the order of (25), i. e. to 2n. By (27) we see that under 
our suppositions the differential equation (25) has 2n linearly inde- 
pendent simple sources 


Etat 


f(r) (7 =1,2,... n). (28) 
As the application of 4,4?, 43... to (28) leaves these fields unchanged 
(so that all radial functions Rj**(r) generated by a given simple source 
(28) differ for different s by a power of —k* only) we have 2n linearly 
independent multipole fields (3) with the same spherical harmonic. 

If one of the kj vanishes, we must replace the corresponding 
exityrlr by 1/r (cf. § 3). 

If equation (26) has a multiple (say »-fold) non-vanishing root 
æ= k", then we have to use in (27) the fields of simple sources belon- 
ging to the differential equation 


(A+k)’u=0. (29) 


Because 


ikr ikr ikr 
ATHENE GG + ikir) & 


and e*"/r is a solution of the equation Au+ k’u=0, we see that a{j+1)- 
fold application of the operator 4+ k? to 7 ett iy yields zero. The 
most general spherically symmetric solution of (29) is therefore given by 


v—1 v—1 


gle e—ikr 
SA ree Goel Sy 
Ud eam > pyri + 2 qyr 


{=o 
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with arbitrary coefficients p; and qj. All multipole fields belonging to 
the same spherical harmonic Y,n(0,p) are here given by the applica- 
tion of the operator on the right-hand side of (13) to the spherically 
symmetric fields 


ikr 


T $ 
in=1 — E O aar D] (30) 


The application of 4,42... to (30) yields namely only fields which are 
linearly expressible by these fields. There are here therefore 2y linearly 
independent multipole fields having the same spherical harmonic. 
If the multiple root of (26) is zero, we have instead of (29) the 
differential equation 
Au =0, (31) 


which has the general spherically symmetric solution 


2v—2 


f(r) = Yar 


=—1 


with arbitrary coefficients cj. Hence, the fields of the simple sources 
are given in the present case by 


f(r)=r1 (j=—1,0,1,..., 29 —2). (32) 


But the application of A,A?... to (32) gives us no new fields of sphe- 
rical symmetry, so that we have altogether 2” such fields. 

As the differential equation (31) is a generalization of the poten- 
tial equation, we must go cautiously to work while determining the 
number of multipole fields with a given spherical harmonic. Applying 
(14) to ri, we get for this simple source 


Rir) =j(j—2)(j—4)... (G—2(I—1)) 9. (33) 


To obtain the radial function Rj**(r), we apply to (33) s times 
the operator appearing in (16). Applying it once torm we get 
(m —1)(m+1+1)r™-*, so that an s-fold application yields 


(m—l)(m—l—2)...(m—l—2s+ 2) 
*(m+1+1)(m+1—1)...(m+1—28+4 3) rms, 
From (33) we get therefore 


Ri*(r) = (§ +1) j(j—1)... (7 28 +2) 


*(j—28)(j —28 —2)...(j 28 —21 4 2) riis, oe 
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By (34) we see that all Rit (r) vanish when j fulfils either the in- 
equality 
—1<j7<2(s—1) (j odd or even) (35) 


or the inequality 
28 SjK28+ 2—2 (j even). (36) 


(35) means that the number of multipole fields Y;m(0,p) Rit™(r) with 
the same spherical harmonic is always finite and (36) that for even j’s 
the number J is limited. For even j’s we obtain therefore only a finite 
number of multipole fields (3) which are moreover all finite at r—0. 
These fields are evidently a part of the regular solutions of (25) of 
the form Y,,(0,~) Si(r), mentioned in § 3 and, in general, not ob- 
tainable by the method used so far. For odd values of j we obtain 
for 1+ 2s<j fields which are also finite at r=0. But for [+2s>j 
they become infinite at r=0. We easily see that this state of affairs 
is a generalization of the one we met in the case of the potential equa- 
tion, where the fundamental laws where given by ri with j=—1,0. 

Concerning the other special cases of (31), we mention here only 
the case of the biharmonic equation A? u= 0. The spherically symmetric 
solutions are here given, according to (32), by f(r)=7?,7r,1,7—1. All 
finite fields at infinity are deducible from the two spherically sym- 
metric solutions f(r)=r and r~t. But it is sufficient to make use of 
the spherically symmetric field f(r)=r only, because then Af(r)=2/r 
and by (18) we have to use besides f(r) also all A*f(r) to get all radial 
functions Ri*%r). 

It may be remarked that (25) seems to be the most general dif- 
ferential equation of a finite order giving fields fulfilling the require- 
ments (a), (b) and (c) of § 1 and belonging, therefore, to an isotropic 
and homogeneous space and to a linear physical law. But in §1 we 
have defined a field without any reference to any differential equa- 
tion, using only a field of a simple source playing the part of a fun- 
damental law. All of our general considerations are applicable to such 
a general case. pa 

If to a given differential equation belongs more than one ele- 
mentary law, we must use in any case all of these laws to get the 
complete field belonging to the given differential equation. But it 
may be also reasonable to define a „partial“ field using only one 
spherically symmetric solution of the given differential equation. 

The question what fields exist with the properties (a), (b) and (c) 
of § 1 which are solutions of an integral equation has not been investi- 


sated in the present paper. 
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§ 5. Fields finite at r=—0. In the case of the potential equation 
the procedure we have used ee not yield fields of the form Y1m(9,9) 
Si(r) with radial functions Si(r) regular at r=0. In the case of equa- 
tion (25) we have obtained a part only of these fields. But we can 
contrive a method to construct such fields also in the general case 
of simple sources which are not defined by a differential equation of 
the form (25). We must only consider the field of the given simple 
source as expressed by a bilinear formula of the sets of both fields 
Vim(9,¢) Ri(r) and Yim(6,¢~) S1(7). 

To achieve this aim we express first the n-th derivative of f(r) 
with regard to z by multipole fields. For m=0 we get by (13) the 
function P0) Ri(r) expressed by the derivatives of f(r), Af(r), A?f(7r) 
with respect to z. If we apply 4,A?, A’... to this equation and lower 
at the same time the index l to 1—2,1—4,/—6... respectively, we get 
the functions 

A(P2Ri2)=P12Ri2, A (Pr Rra)=PraRa, 


A (Pire Ring) Prs Rie. 


expressed by the z derivatives of Af(r), A?f(7), A?f(7r)... Multiplying these 
relations by suitable constants and summing up we get finally 


Jr [n/2] vps (n—-==) az 
= =y! i i 2 Rn—s(r)Pn—s(0) (37) z 
en > (2(n—2s))! 2s—3 24s g! ; 


Denoting by o the distance of a point P, situated on the z-axis 
from the origin O of Cartesian coordinates, by r>o the distance from O 
to the point P(x,y,z), and by R the distance from P, to P, we have, 
as R= ?+ y?+ (z—o) and r?=g+y?+ e, the Taylor series 


OS Eh (38) 


Since 


221 T( wep a 


we get from (37) and (38), putting n=l+ 2s and using J and s as sum- 
mation indices, the relation 


}(R)= Sc =i pees Pi MO Va Srey ga = Oe 


War s+) \2 


we intended to obtain. 
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First we show in special cases that (40) yields well-known for- 
mulae. In the case of the potential equation (f(R)=1/R) Rir) is given 
by (19) whereas all other Rit(r) vanish. By (39) and (40) we get 
therefore the generating function of the Legendre’s polynomials 


1 SPO) 
SDN rT e (41) 


In the case of the wave equation Au+ u=0 we have f(R)=i e IR 
and according to (22) and (23) 


Bir) = (—1yits FA) 
r 


Because 


co 


a Me 1 o\H2s+1 
vie) = Vr > ( 1) s!T(l+s+8) (5) ; (42) 


s=0 


we get the well known relation due to E. Heine 


eR 4-1 Clr) yilo) o- 
Nr E = EMA (*) 


E= 


The Cartesian coordinates %,y,z used for the establishment of 
relation (40) do not appear any more in this formula. 6 denotes here 


> > > > 
simply the angle between the vectors OP»>=o and OP=r. Denoting 


> > 
now by 0, 0o, Po and 7, 6,, pı the spherical polar coordinates of @ and r 
in an arbitrarily fixed system of reference, we have by the addition 
theorem 
+ i! 
a ' * 
P{0)= 2 Tmt =m Yim(9o,Po) Yim 91,71), (44) 


m=— 


where the asterisks denote the conjugate complex., values. If we in- 
troduce (44) into (41) and (43), these formulae express Green’s function 
for infinite space for the potential and the wave equations by means 
of a bilinear combination of two sets of solutions of these equations. 
In the case of (41) these solutions (without taking into account the 
normalizing constants) are given by Yim(9,%)e’ and Yim(6,,9,)/7", 
and in the case of (43)¢by Yim(Oo,%) ye e@)/e and Yim(%,91)6(7)/r. 
Without any reference to the respective differential equations, 
we get the finite solutions Yim(9%,%)e! and Yim(91,9,)yi0)/e as fol- 
lows. We multiply (41) or (43) by Yr,m(61,91) and integrate the pro- 
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duct over the unitary sphere. This procedure is equivalent to sum- 
ming up the fields f(R) with differently situated source points P, and 
therefore must also lead, like derivation, to ,,fields” in the sense of 
the definition given in § 1. In the same way we get from (40) 


eS RH s : 
TAs or Po) Dyer ie (45) 


where we have written again l,m instead of V’, m”. 

All expressions (45) with different l,m are linearly independent 
because they contain different spherical harmonics. But they do not 
depend only on ọ, 6, and pọ but also on r. To obtain fields of the form 
Vim( 95) %o) Si*’%(o) (in the sense of §1) depending only on @, 4, Po 
let us suppose that all the R} (r) with a given l are expressible by 
a finite number of them, say n. Expression (45) splits then into n 
groups each of them multiplied by one of the linearly independent 
R (r). As (45) represents (as a function of o, 0o, Po) a field for every 
special value of r, each of these n groups may be considered as re- 
presenting`a field. In such a way, we get n fields YVim(0o, Po) St Co), 
each of which is expressed by a power series and is finite for o= Q. 
The numbers of the fields Yim(0,91) Ri (r) and Yym(9,%) S He ) 


belonging to,the same indices l,m are generally the same. Because ~ 


the fields Yzm(0,~) S}(r) are everywhere finite in the finite domain we 
can imagine that, in agreement with the above considerations, they 
are generated by elementary sources situated at infinity> 

The method used in this section rests in reality upon the fact 
that it is possible to define by an elementary law er a simple source 
a field w outside a system of simple sources distributed over a space D 
by integrating over this domain the product of the gemon EG law 
f(E) and the strength g of these sources: 


u= [f(R)gdr. } os he} 
D 


For f(R) we may use here expression (40), supposing that the point O 
is situated either inside or outside of the spatial domain D. If O lies 
within D we have to integrate over the variables ọ, 0p, pọ and obtain . 
for u an expansion in the multipole fields Yym(6,, pı) Nitto A Wee But 
if O lies outside.D, we have to integrate over T, 01, pı and obtain u 
(46), as a function of 0; 9, Po, possibly as an expansion in a series of 
the functions Deen E EN (o) if they exist. 
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In this connection we may ask what fields 
a [RO(R) go at (47) 
D 


we obtain if we use in (46) instead of the elementary law f(r) any sphe- 
rically symmetric field 


A* f(r) =Ro(r). 


Physically we have to interpret these fields outside the domain D as 
the field of a spatial distribution of spherically symmetric multipole 
sources (cf. Appendix) generated by systems of multiple sources of 
order 2%. By Green’s theorem 


2s 2s—2 
[eaa RE) ae S Feat Sj z do 


On on 
D 


and we see that 
2s—2 
us = eh Ro? (R) A go’dr — f (z; sA a e gp 2R E) do. (48) 


It is therefore possible to reduce the integrals (47) finally to the vo- 
lume integrals (48) with the elementary law f(R) and a sum of sur- 
face integrals with Ro” (R), oie ... We may of course also define 
fields i 
| ui = f Yum(0,p) RR) gi dr 


of spatial distributions of multipole sources with a directional sym- 
metry given by the spherical harmonic Y7,,(6,y) and reduce them in 
a similar manner by Green’s theorem. 

But, we must remark that it depends upon the order of the sin- 
` gularities of f(r) or Ri (r) at r= 0, whether the fields u (equation (46)), 
or uit respectively are finite within the domain D of the source dis- 
tribution in question. In the potential theory, for example, already. 
quadrupole sources give rise to infinite fields in the domain D. We 
have therefore to exclude, as the case may be, some elementary laws 
if from the physical point of view we must suppose that the field wu 
or one of the fields uitis: finite in the domain D. 


§ 6. Fields of simple sources which are not solutions 
of the differential equation (25). In § 3 and § 4 we have studied 
fields of simple sources which were solutions of the differential equa- 
tion of finite order (25). We can therefore decide now which spherically 
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symmetric fields are solutions of (25) and construct such spherically 
symmetric fields which do not solve this differential equation. 

As we shall see later on in § 8, we are interested in this case in 
the question how many fields 


Yim( 9,9) Rit r) (49) 


exist with the same spherical harmonic Y,,,(9,y). We obtain all the 
radial functions R}**(r) belonging to a given fundamental law f(r) by 
applying the operator in (13) to the spherically symmetric fields 


1) Af(r), ASe (50) 


The number of linearly independent fields (49) with given indi- 
ces l,m is surely finite if 
(A) between the spherically symmetric fields (50) there exists 


a linear dependence of the form (J; a, 4°) f(r) = 
s=0 


(B) the number of the spherically symmetric fields (50) is finite. 

In both cases we have to do with that partial field of (25) which 
belongs to the given elementary law f(r). In case (B) we have 
A’f(r)=0, which is a special form of (25). 

A fundamental law f(r) which is not a solution of the differential 
equation (25) cannot therefore fulfil condition (A) or (B). Consequently, 
we expect that for an arbitrarily given fundamental law there exists 
generally an infinite number of multipole fields (49) with a given sphe- 
rical harmonic. 

Finally, we quote examples of fundamental laws not solving the 
differential equation (25). From the considerations in § 4 we see that 
ri or rie*/r are such fundamental laws if j is an integer smaller than 
—1 or an arbitrary (complex) number not equal to an integer. Com- 
bining such fields with each other or with spherically symmetric so- 
lutions of (25), we can build up a great variety of such elementary laws. 


§ 7. Other applications of formulae (13) and (18). Our re- 
sults given in § 2 are applicable also to time-dependent scalar fields. 
As an example we consider the field which is a solution of the time- 
dependent wave equation 

Au= t a 
E JE 6D) 


The field of a simple source is here represented by 
F(r—ct) 
r 


G(r—ct) 


f(r,t) = e =o | (52) 
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As together with (52) all the derivatives of f(r,t) with respect to LYZ 
are also solutions of the wave equation (51), we can build up solutions 
of (51) having the form of multipole fields (3). We may therefore 
apply (13) to (52), and so obtain solutions of (51) of the form 

z 1 d\'/F(r—ct) G(r—ct) 

Vim(Gyp)r'(= 5) AEA 5 e), (53) 
This result is known (Forsyth 1912, p. 496 and 900) but our deriva- 
tion is considerably shorter than the usual one. Application of 


È a? 
In the case of the differential equation of heat conduction we 
may also apply formula (13) to the spherically symmetric solution 
corresponding to an instantaneous point-source of heat in order to 
obtain fields of instantaneous multipole heat-sources. 
Another application of formula (13) and formulae (18) is the 
following: Let for example G(S,P) be a Green’s function (with the 
.source point S and the field point P) for the potential or the wave 
equation, belonging to the boundary condition u=0 or S40. Then, 
applying the operator on the right-hand-side of (13) to the coordinates 
of the source point S in G(S,P), 4sG(S8,P), AsG(G@,S)... (As is the 
Laplacian with respect: to the source coordinates), we get fields ful- 
filling the given boundary conditions and having singularities in S 
given by the multipole fields Y;m(0,p) Ri (r) (s =0,1,2...). 
Thus we can for example, starting from a Green’s function 

construct a velocity potential field which is a solution of the wave 
equation Au+k’u=0 and corresponds to an arbitrary acoustical source 
(a loudspeaker) given by a series of multipole fields. 


§ 8. Fields defined by elementary laws. To solve the boun- 
dary value problem for the potential or the wave equation in the case 
of a sphere we may use the solutions of these equations obtained from 
the field of a simple source or by the procedure described in § 5. For 
the space inside the sphere we have to use the solutions Yzm(9,o) @! 
and Yim9,%)yie)/o regular at 9 =0 and for the space outside it 
the solutions vanishing at infinity and given by Yj,n(,,9,)/r'+! and 
Yim(91,91) €(7)/r. For a space between two concentric spheres both 
sets of solutions must be used. In all these cases we may choose as 
boundary conditions on the surface of the sphere either the, values 
of the fields or its derivatives in the direction of the normal. 

The application of explicit solutions guarantees that all peculari- 


ties of the boundary problems in question are automatically taken 
12* 


10" \? : : 5 ee : f 
* = ( ) to (53) yields the time derivatives of this expression only. 
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into account. In the case of the potential equation the potential func- 
tion remains undetermined to an additive constant if the derivative 
of the field in the direction of the normal is given on the surface of 
the sphere. If in the case of the boundary value problem of the wave 
equation for the interior of a sphere the radius of the sphere happens 
to correspond to an eigenvibration, then the amplitude of this eigen- 
vibration remains undetermined. We may also require that on the 
surface of the sphere a higher, say the n-th, derivative of the field 
taken in the direction of the normal has a given value. For such boun- 
dary conditions, however, some of the amplitudes (in the case of the 
potential equation the a n—1) may remain undetermined. 

The solutions of the biharmonic equation Au=0 which are 
regular at infinity are e E as we have seen in § 4, from the 
field of the simple source f(r)=r. As Af(r)=2/r+0 and A*f(r)= 
only the functions Ri(r) and RP (r ) do not vanish, whereas all other 
Ri(r) are zero. Two radial fanctions Ri(r) and Rj*(r) belong to every 
spherical harmonic Y,(6,~), so that we have to assign now two boun- 
dary conditions:on the surface of the sphere, e. g., for the field and 
its first derivative in the direction of the normal. 

In the case given by the solutions of the differential equation (25), 
there belong to a given spherical harmonic Y,,,(6,g) generally n linearly 
independent radial functions which are finite at r=0. We can there- 
fore fulfil in the case of the interior boundary value problem n boun- 
dary conditions on the surface of the sphere, e. g., we can assign the 
value of the field and its n—1 derivatives in the direction of the normal. 

The complete system of the solutions of a differential equation 
in spherical polar coordinates being given, we can expand any other 
solution into a series of these special solutions. By the help of such 
expansions we can solve, in principle, every boundary value problem 
also for a space of non-spherical shape. 

As we have shown in this paper it is possible, however, to build 
up fields of the form Y,,(6,p)Ri (r) and Y,m(0, p) S(r) starting 
from the fields of onẹ or several given simple sources (or elementary 
laws, which may be given empirically). We may therefore ask if there 
is any possibility to formulate and solve boundary value problems in the 
domain of a field which is not given by a differential equation of the form 
(25). If the space of the problem is-a sphere, we may use for the solu- 
tion of the boundary value problem the multipole fields (3) and the 
finite fields Y,m(0,p) S} (r) of § 5. If the number of the fields be- 
longing to a given spherical harmonic is finite, say n, and n is inde- 
pendent of the order of the spherical harmonic, tien on the surface 
of the sphere we may assign arbitrary values to tHe field and to n—1 
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of its derivatives in the direction of the normal. It seems therefore 
that also in the case of a space of any shape it is reasonable to put 
boundary value problems only if the number of fields belonging to 
one spherical harmonic is finite and independent of its order l. But 
so far no example has been found for such an elementary law not being 
a solution of the differential equation (25) in concert with our con- 
siderations in § 6. It seems therefore that we must consider as open ~ 
the question whether there exist fundamental laws not being solu- 
tions of (25) but giving fields in the domain of which we can put boun- 
dary value problems. 

It may be possible, however, to formulate and solve eigenvalue 
_ problems for fields given by a fundamental law only, e. g:, for a sphere, 
if the field of ‘the simple source contains an arbitrary parameter k, 
as in the case of f(r) = g(r) (e#—e—*)/r. The function g(r) must of 
course include also negative powers of r (ef. (30)), otherwise we 
would have to do with a solution of the differential equation (25). 
For the boundary condition to which we may subject the field it 
is sufficient to take u=0 or du/dn=0. 

In the present paper we have considered scalar fields only, but 
vector, tensor, and spinor fields may be treated in a similar manner. 

It was not possible to answer here satisfactorily the question 
which problems can be solved by using fields defined by fundamental 
laws only. To do this would require a thorough mathematical investi- 
gation. But it seems desirable to follow up this question out of regard 
to the possibility that the differential equations which form so far 
the basis of all our considerations in theoretical physics may turn 
out to be not sufficient for representing the physical reality. In such 
a case one could try to use the more general fields defined by funda- 
mental laws. But even if it should prove to be certain that we do not 
need these more general fields in physics, their study shows us clearly 
the limitations imposed on physical fields by the differential equa- 
tion (25). We may also remark that by using elementary laws we get 
the possibility of estimating the effect on the field caused by smaller 
or greater changes of these laws. 

There is no need to emphasize that we can also find elementary . 
laws giving fields fulfilling the requirements of relativity theory. 
A simple example is the elementary law (52) yielding fields which are 
solutions of the time-depending wave equation (51). Finally we may 
remark that boundary*value problems are very important in classical 
but not in quantum physics, and the possibility or impossibility to 
solve such a problem in the domain of a given field does not generally 
decide about the applicability or not of this field in quantum theory. 
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Appendix 


about the systems of simple sources producing fields 
of multiple sources and multipole fields 


. v 
Maxwell (1892) defines the field of a multiple source of order 2 by 
at f(r) ae 
Ne 08; 08g... Oy’ 


where A 3 


28; a teas Loum eee 


denotes a derivative along the direction given by the cosines aj, pi Yi. 

The field (A1) is expressible by a linear superposition of the spe- 
cial fields (1) containing 21’'+1 constants, 21’ corresponding to the V 
directions of differentiation and the constant A determining the am- 
plitude of the field. A field given by an arbitrary linear superposition 
of the fields (1) is generally representable by 


f(r) Pert mee ore ee ae 


Tee O84 ISg...08y 0840983... 0Sv_2 Os Oso ... O08" —4 


se. (AB) 
because the numbers of arbitrary constants on both sides of this equa- 
tion are the same (år). 

Assuming that the simple sources belonging to the systems of 
simple sources of the different fields (1) do not generally annihilate 
each other, we can state that the field we get by a linear superposition 
of the Ay fields (1) of multiple sources of order 2” is generated by a sy- 
stem consisting of Ay 2” simple sources. On the other hand, we may 
interpret the field on the right-hand side of (A2) as due to a system 


of [2] +2) 2” simple sources. The number of simple sources we can 


imagine to produce a given field obtainable by a superposition of the 
fields (1) with a given J’ depends therefore on the representation we 
use and is not determined uniquely by the field itself. 
A simple example of this statement showing its mechanism is 
f(r) | 2 f(r) ne) of(r) 

Alora, Rae Hae |e cee 
where the same field can be interpreted either as due to three double 
sources or to one double source only. Here we can clearly see. that 
the simple sources of the three double sources on the left-hand side 
of equation (A3) do not annihilate each other, but the vector addition 


of their dipole moments gives the dipole moment of the double source 
on the right-hand side of this equation. 
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Which system of simple sources produces a given field can be 
decided on the basis of physical arguments only. 
The fields 


f(r) =A'f(r) 


which present themselves in (A2) are in spite of their spherical sym- 
metry not produced by simple sources, but by a superposition of fields 
of multiple sources of order 2”, 

It follows from the above considerations that the number of 
simple sources constituting the source system of a multipole source of 
order 2! is generally not given by 2 but determined by the order 2” 
of the fields of multiple sources of which the multipole field is built up. 

If a field is a solution of the potential or the wave equation or 

in general of the differential equation (25), there appears a new in- 
determinateness in the system of simple sources assignable to a given 
field. Putting field (1) in one of these equations, we obtain a vanishing 
linear superposition of these fields. As their systems of simple sources 
do not generally annihilate each other, we obtain in such a way source 
systems producing no fields. It is possible to add such a source system 
with a vanishing field to the source system of any field without chan- 
ging it. It is therefore generally impossible to assign uniquely a system 
of simple sources to a field which is a solution of the differential equa- 
tion (25) if there is no additional physical information available about 
the source system of the field in question. On the other hand, it is 
possible, in principle, to ascribe also to every regular point of the 
field such a source system, and to the whole field an „infinite sea” 
of such source systems. 
If we do not make any explicit use of this possibility and assume 
that the source system of a field is given by the derivatives which 
have been employed for obtaining it from the elementary law, we 
have to make the following remarks concerning the source systems 
of the fields which are solutions of the differential equation (25). 

In the case of the potential equation, the source system we can 
ascribe to the multipole field (20) with given indices l, m depends 
upon the choice of a, b, c. i i 

In the case of the wave equation Au + kèu=0, we have 


u=(— 5) 4u so that we can interpret for example the spherically 


symmetric field et#r/r by assuming that it is caused either by a simple 
source or by a system of multiple sources of any order ek 

In the case where the 2n fields of simple sources of the differential 
equation (25) are all different and given by (28), we may make the 
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same remarks as in the case of the potential or the wave equation. 
Of a particular interest, however, are the solutions of (25) which are 
at the same time also solutions of equation (29) or (31). In these cases 
we may interpret all the spherically symmetric solutions (30) and (32) 
as generated by simple sources. But it is also possible to assume that 
only two of all these spherically symmetric fields are generated by 
simple sources, whereas all the other fields (30) and (32) are due to 
higher multiple sources. In the case of equations (29) and (31) we can 
assume that the fields of simple sources are given by 1-1! et#r/r and 
p22, 2-3 respectively. All other spherically symmetric fields (30) 
and (32) are namely obtainable by linear combinations of spherically 
symmetric fields which we get by repeatedly applying the operator 4 
to the above mentioned fields of simple sources. 

If we regard, as in quantum theory, the field as generated by 
some elementary particles, we have generally to ascribe different fields 
of simple sources to different particles. Their minimum number is 
therefore given by the number of different roots of equation (26). 
In the case of a simple root it is namely possible to assume that one 
of the spherically symmetric solutions e+j'/r (with a positive or nega- 
tive sign) is due to one elementary particle and the other one is pro- 
duced by particles which are at infinity (cf. § 5). In the case of a non- 
vanishing multiple root of (26) we can similarly ascribe one of the 
solutions 7” et*r/r to one particle and imagine that the other spheri- 
cally symmetric solutions (30) are due to multiple sources or to par- 
ticles situated at infinity. If the multiple root becomes zero, it is suf- 
ficient to regard the solution r?”—! as the elementary law, the other 
solutions (32) with odd indices as due to multiple sources, and those. 
with even indices as belonging to solutions which are finite at r=0 
and generated by simple sources at: infinity. 
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A NOTE ON PERTURBATION THEORY 


By Jan RZEWUSKI, Physical Institute, Nicholas Copernicus 
University, Torun. 


(received August 8, 1951) 


In this paper the identity of the results of the new and the old perturbation 
theory is shown by a direct transition. Feynman’s method is presented in a simple 
form based on integral equations which is evidently invariant not only in the results 
but throughout the derivation. The problem of scattering on bound states is treated 
by means of the modern techniques. 


Introduction 


At the present stage of perturbation theory one has to consider 
two methods. The old one described e. g. in Heitler’s book (1947) and 
‘the new one developed by Schwinger (1948, 1949), Feynman (1949) 
and Dyson (1949). The difference between these two methods les 
apparently in their invariance properties. The old method starts with 
the non-invariant Schrédinger equation and uses a formalism of field 
quantization based on non-invariant commutation relations. The new 
method starts either with an invariant generalization by Tomonaga 
(1946, 1947) of the Schrödinger equation or with a Langrangian, and 
uses either a formalism of field quantization based on invariant com- 
mutation relations or a semi- -classical treatment based on Stueckel- 
berg’s causal A, function. 

Feynman (1949) has shown that one may derive the new results 
also from the non-invariant Schrödinger equation, at least in the case 
of a given external potential. As the two methods start with the same 
Lagrangian and both use an expansion in powers of e, their equi- 
valence seems to be proven beyond any doubt. The results, however, 
seem at first sight to be incomparable since the formulae are con- 
structed by means of different quantities. In the old method the ele- 
ments of the perturbation formulae are the transition amplitudes and 
the energy differences, in the new method the formulae are construc- 
ted out of various types of invariant A-functions. 
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It is the subject of this paper to show, that both these seemingly 
different forms of perturbation theory are identical. The argument 
consists in a direct transition from the results of the new perturbation 
theory to the results of the old one and vice versa. In other words, 
it is shown that one could obtain the covariant formulae of Feyn- 
man and Dyson without using their covariant formulation, but simply 
by carrying the old theory just one step further. 

In the first section the prescriptions of the new perturbation 
theory are shortly deduced by the use of Feynman’s semi-classical 
treatment but without its non-invariant aspects. In the second section 
the transition to the old form is explained on three examples. Finally, 
the results of the calculation are discussed. 


§ 1. Covariant perturbation formulae 


To introduce the covariant perturbation formulae we choose 
a formalism based on Feynman’s semi-classical treatment. To avoid 
its non-invariant aspects, however, we consider integral equations 
rather then differential equations and introduce a covariant definition 
of the transition amplitude. 

We start from a Lagrangian describing the electron field in inter- 
action with the electromagnetic field 


L=y(tV—m)yt+ tl wl w—epAay, (1) 
where 
eA a 
Pw Fo, Oily’ 


A = ypÅpu= YÁ — yA — yz Áz — yÁ, 

k Oat kee ð ð ð 
N. — Ye, aa aa igy aga 
En Rs T 2 
OL, EN da,’ da,’ ey 


For particulars of this notation see Feynman’s (1949) paper. The dif- 
ferential equations following from (1) are 


> < 
(tV—m)y=eAy, pee (3) 


OAp=evypy, 


where the arrow indicates the direction of differentiation. We re- 
place (3) by the corresponding integral equations: 
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pL) =p(1)—ie | K(12).A(2) y(2) dary, 
P) =p (1) +40 f dayp(2 TNT 
BEDEN T e [PDP y y2) da, 


dx, = (dX; di,dx3dx,4)x=x,, 


(4) 


It is easily seen that the functions y and A, satisfying (4) are solu- 
tions of (3) if 


(iV —m), K(12) = iô(12) (5) 
K°(21) = BK(12)*8 
(14°(12) = 6(12) (8) 


It follows from (5) that K° must satisfy 


zs < 

K°(21)(tV—m),=76(21) (5°) 
The suffix „1” denotes that the corresponding operator acts on the 
variable æ. The functions y®, y?, A®° are solutions of the corresponding 
unperturbed equations 
. > < 

(V—m)yP=0, y(iV+m)=0, 04%=0 (T7) 

The third equation (4) has already the desired form since it expresses 


A,, in terms of A‘, @, y and y. The first two equations contain y in 


the integrand and outside the integral. To solve them for y one has 


to expand in powers of the small coupling constant e 
yo. 
yO(1) = —te a K(12).A(2) y°(2) da, ee 
p(1) = (—te af ‘[®(12) 12) A(2) K(23).A(3)p(3)da,dars 


SEtelete =.s) eule. te) ioe 6 ey vm 16 ee pe en) al. té~ ,0) Oo of ona ce A a er ie 


Summation of (8) yields 
= XJ y(1) = y(1) te | K(12) A(2) (2) dary (9) 
n=0 


with 
K(12)=K(12)—ie fE K(13) A(3)K(32) datg+ ... (10) 


Consider now for a moment the field A as given. In this case, 
we may express the scattering amplitude for one particle by means 
of the covariant expression 


R= [pol (iV—m)y, du 
co fgeaptan—ietf f7) K(12).A(2) p%(2)da, dæ, 


(11) 
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where q and p are the initial and final momenta of the particle. K? is 
well defined by equation (5) and the demand that it describes pair 
production and annihilation. The connection of (11) with the con- 
ventional transition amplitude S is 


y i ox 3 ESMI } aaO ) 12 
S= f' pty e= [yiye —iR= (2x p—g)ðy— iR (12) 


Here i, j numerate the four solutions of the Dirac equation corre- 
sponding to the same momentum. The ôy factor expresses the con- 
servation of spin and of the energy sign. 

To make the transition to the old theory we shall need a rela- 
tion between K? and the eigenfunctions. For discrete eigenvalues this 
relation is given by Feynman (1949). The generalization for continuous 
eigenstates is immediate. Indeed, the solution of the first equation (5) 
may be written 

4 e—tkx12 


A2) = Cat! eam 


(13) 


To account for hole theory the integration path must pass around 


the negative pole pee pee m? in the lower and around the posi- 


tive pole ka= + Vie m* in the upper half of the complex k,-plane. 
Carrying out this integration, we get 


1 sp Etm 


K?(12)= —iRkX 45 
where the upper sign corresponds to 4, >t, eee ake m?, and the 
lower sign to t <t», k= Vtm. Now, it is easily seen that 
2 3 
> for k,>0 
k+m\ _ Jit |, orno 
Diels me eal eee ee (15) 
D for k>0 


i=3. 


The upper indices 1 and 2 correspond to ky>0, 3 and 4 to k<0. In- 
troducing (15) into me one gets 


Pier aan Sa ‘wOlk,1)ph(k,2) for t >t, 


rials a : (16) 


-Sa p f ee ST le oraa 


z of 
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Now, we may go on to consider à special kind of transition in 
which the final state is the same as the initial one, p=q. The first 
term in (11) must be omitted, otherwise p=q would imply no scat- 
tering at all and therefore A= const. The second term can be sum- 
med for t <t, over the positive energy states, for t, >t, over the ne- 
gative energy states. Using (16) one gets then i 

2 


+> | 
“i Ete ff p°(p8) AL) K (12) (2) "p(p2) da, da, 
(17) 


=—ieSp | f K(21)A(1)K(12).A(2) dey, dirg.. 


This is the probability amplitude for the creation and annihilation 
of the same pair of particles. We call (11) an open and (17) a closed 
cycle. 

Interaction between particles by means of the electromagnetic 
field is introduced formally as interaction between a number of such 
cycles. One has to replace each pair of A(a;)....A(a,) by a factor - 


1 
Yu: Yu GM Uide), (18) 

where @(12) is a solution of 
(,G9(12) = 6(12). (19) 


Each possibility of combining the A’s into pairs gives rise to an extra 
term in the transition amplitude. The factor (18) corresponds to the 
emission and reabsorption of the same virtual quantum. Emission or 
absorption of a real quantum is described by a single factor A sui- 
tably normalized. It may be noted that a cycle having not less then 
two factors A may also interact with itself. A closed cycle with an 
even number of A’s interacting only with itself is unobservable and 
has, therefore, no physical significance. Examples will follow in the 
next paragraph. 


§ 2. Transition to the old perturbation formulae 


Having -thus introduced covariant prescriptions for calculating 
transition amplitudes, we have to show their connection with the old 
formulae. It is sufficient: to carry out the calculations for some typical 
examples, since for more complicated problems they run in exactly 
the same: way. We choose the electron self-energy problem, the Comp- 
ton effect and the problem of scattering on bound states. 
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(a) Electron self-energy. 

We treat the problem in first approximation. This means one 
open cycle with two factors A interacting with itself. According to 
(11) and (18) one has 


Rear=— f [PPI yu K2) yuy?(q2)G(12)dar,da,. (20) 


This formula is evidently covariant. Now we have to spoil this co- 
variance in order to make the transition to the old perturbation theory. 
For K?(12) we insert (16) and for G°(12) a formula analogous to (14): 


3 3 
G°(12) a ee d* k 


a —ikry. € 
(27) J 2k 4 i eee 


> 
where the — and the + sign correspond to t >t,, ka= + |k| and t< tə, 


> 
k,=—|k| respectively. This form of G° is necessitated by the physical 
interpretation. With (16) and (21) we get for the self-energy 


2 
Dyusfobint, Sipe 0 


pense Pr etry) 
self (2218 Ly ALa if 4 


oe Lerten ge 
s=3 
= 1 
X YAPI) yup (r1): P(r) Yad p"(q2) g ee. 
ag 
Inserting plane waves for the y’s 
Fy q 2) = 4g) e—a% (23) 
and separating the time variable, one gets 
2 
> t>tk,>0 
Ramie | @rark | = Nema (cg) 
D> tis hky<0 
3=3 á (2a 
> > J Te >w abs 
x Çi, plyuls,7><s,7 Yali g>: ok, e'h Orki tka] dt, dt, 
x 4 ? 
where ny 


> > 
Li, ply ul 87> = PAP) (Yap Yr). (25) 
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It is convenient now to change the time variables: 


te =t, arose 


Bamia faf Se fe +(e ie et (26) 


ge “te | See y ; 
x <i, plyus, p —k> <8, —hlyal i, p> = x f eae PAR 
hee 


Detailed calculation would show that terms with ij vanish. Carrying 
out the integration over t and t, we may thus finally write 


ent Rig ie 
dk Sik, plyals, p — k>} ÈKi Pirs pi> 


is ee es r e A E 
[k] pa—lb|—Vip— hy? m? pat [kM (pk me) 


This is exactly the self-energy formula as calculated by the old per- 
turbation theory. 


o 


(b) Compton scattering. 

The formula for Compton scattering is given, in first approxima- 
tion, by the second term in (11) with K replaced by K° and A pro- 
perly normalized 


Vat —ie f fy (p1).A(1)K(12).A (2) y(q2) da, dz, (28) 


It may be easily shown that A must have the form 


A(z) = re (ee gikx) cared E cos (kx) (29) 
V2k, Ky 

if we demand that it represents one photon in a box of unit volume. 

It is sufficient to consider only one term of (29). The contribution 

of the other is then obtained by a suitable change of sign in the kps. 

Putting also (16) for K® one gets 


O kle® 
R,=(2 z} (p+ k—1—q) == BOSE „ple |s, p+ ED <8, Bt kje ind (30) 


ka Vikas Pat ky— 
with 
l E,=Er =+ Vp +k m, (30°) 


B= By = — Vpt k+ m. 
This is the same formula as that derived by means of the old pertur- 
bation theory. 
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(c) Bound states. 

In the problem of bound states the method of § 1 is applicable. 
The only change is that solutions y® of the unperturbed Dirac equa- 
tion are now replaced by the solutions of the equation 


(iV—m+ V) yp=0, (31) 


where V, is any potential allowing a discrete spectrum of energy eigen- 
values. The formula (11) for the transition amplitude is valid with y® 
replaced by Yn and K? given by an equation analogous to (16): 


2 Pap) for t, >t, 
(12) = (32) 
0 m SO ila. 


One puts 0 for t,<#, to avoid creation of antiparticles. When V is 
stationary the y, have the form 


> 
Yn(@) = Yr(©) en (33) 
and (32) becomes 
> > 
D Pnl) Pnl Ta) e— iEn —h) for t >t 
man =| (34) 
0 TA TAN 


The factor A occuring in (11) represents absorption or emission of 
real or virtual photons respectively. In the case of real photons we 
have to replace it by (29). We take again only one of the summands 
n (29). Considering in (11) only the two terms proportional to e and e 
one gets with the help of (29), (33) and (34) 


R= a YmAypnda—ie 5f i Pm(1)A(1) Spl )p2)-A (2) pn(2) dar, dezt... (35) 


t >t 


Separating here the time from the space integration one gets in the 
same way as before 


R=22¢65(Em—hy—En j <mlect##|n> 
V2k, sF (36) 
+ 26 5(Em+ ky qa —En) AET PAG 
| 4kaqa(Em—qa— Er) 
where ; 
Ss joe ere’ a > 
<m eze ny = [Pm a) epii ya) Bar (37) 


The first term gives the probability amplitude for the absorption of 
a quantum k accompanied by a transition of the bound system from 
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the state n to m. The corresponding emission amplitude is obtained, 
apart from a weight factor, by changing k to —k (second term in (29)). 
The second term in (36) gives absorption of q and emission of k. Other 
processes of second order in e could again be obtained by taking all 
possible combinations of sign in the q and k 4-vectors. One would 
thus obtain the amplitudes for emission of q and absorption of k, 
emission of both q and k, and finally absorption of q and k. 

On these examples one can see that it is possible in any case 
to go over from the new to the old formulae by the use of (16) or (34), 
(21) and by separating time and space integrations. Of course, the 
reverse transition from the old to the new formulae is also always 
‘possible by going the same way in reversed direction. Thus, this is 
a method to obtain the new covariant formulae by means of the old 
perturbation theory. 


Conclusion 


The calculations of the preceding paragraph show directly how 
the non-invariant conceptions enter into the theory. It is the space- 
time integration carried out in the transition from (22) to (27). The 
integration over space gives 6-functions depending on the differences 
of momentum vectors in the successive states. These 6-functions intro- 
duce conservation of momentum. The energy is conserved only in 
transitions between two real states. For transitions between two vir- 
tual or one real and one virtual states there is no conservation of 
energy. This is caused by the discontinuity of the kernels K? and G° 
at a hyperplane t= const. The conservation of momentum without 
conservation of energy is a non-invariant feature since it can be valid 
in one frame of reference only. In the new formulae this feature is 
removed by putting energy and momentum on the same footing: 
Both, momentum and energy are conserved in transitions to virtual 
states. The corresponding four-vectors, however, are not restricted 
by the constancy of their squares. This restriction is replaced in 
the new theory by the invariant property of the kernels that’ their 
only contribution stems from the residua at p2—m?= 0 or k?=0, where 
p and k are the momenta of the intermediate electron oF. photon re- 
spectively. È 

There is one further question which we should like to discuss 
here. It is the fact that the new formulations of the perturbation theory 
by Tomonaga, Schwinger, Feynman, Dyson and recently by Yang 
and Feldman (1950).do not bring anything new from the point of 
view of the invariance properties. As already stated the covariance 
was contained in the old formulae and it was only necessary to make 
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one step to bring it to light. Some of the new techniques, however, 
have other merits which distinguish them from the others, as well — 
as from the old theory (we do not wish to discuss here the ideas of 
mass and charge renormalization, which stand quite apart from the 
problem under discussion). Such are the techniques of Feynman and — 
Yang—Feldman. Indeed they permit quantization of fields (Rzewuski — 
1950, Rayski 1950) which could not be quantized by means of — 
the old theory and the formulations of Tomonaga, Schwinger and — 
Dyson. An example is the field of an extended electron. Both ~ 
techniques achieve these ends by introducing a method of interpre- — 
ting quantum statistically a formalism in Lagrangian form without 
the help of the Hamiltonian formulation. 
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MEASURING ENERGY EMITTED BY ELECTRICALLY 
ACTIVATED NITROGEN * 


By H. CYGAN, Physical Laboratory II of Wroclaw University, 
Wroclaw 


(received June 27, 1951) 


Lord Rayleigh’s work on the measuring of energy emitted by electrically 
activated nitrogen was repeated, and the causes of the very high values of energy 
he obtained were studied. One comes to the conclusion that the results obtained 
by Lord Rayleigh may be easily explained. 


1. Introduction 


This work was undertaken to verify Lord Rayleigh’s results (1947). 
Very high values of energy were obtained by him from a certain vos 
lume of electrically activated gas (H,, O, Na). The following result- 
were established by him for nitrogen: 


pressure in the discharge lamp (mm Hg) 0,26 0,124 0,052 0,026 0,015 0,011 
Energy in eV for every particle in the lamp 2,16 4,37 21,1 70,6 147 233 


The method of measuring and the construction of the lamp remain 
essentially unchanged. 

A certain volume of gas in a cylindrical discharge tube was acti- 
vated by an electrodeless discharge by means of impulses of a high 
frequency current, each lasting several microseconds, passing through 
a coil wound around the discharge tube. In the period between two 
consecutive impulses (0,5 sec) nitrogen diffused to a measuring device 
placed in the axis of the discharge lamp. and there lost its energy in 
form of heat. The measuring device consisted of a band of platinum 
foil whose temperature, was controlled by measuring its resistance in 
a Wheatstone bridge system. 


* Reported in the XIIIth Congress of Polish Physicists in Cracow (1950). 
13* 
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A thinwalled glass tube could be placed on the foil to prevent 
the activated nitrogen from coming into contact with it. Measure- 
ments were made with the foil screened and free. Energy from the 
gas when the foil was screened was replaced by a stronger current 
through the bridge, one branch of which was the foil. The resistance 
of the foil during each measurement was maintained constant. 


2. Description of the apparatus 


The apparatus was composed of four parts: 


(1) The measuring lamp with a vacuum set. 


(2) Electric apparatus to activate nitrogen by impulse of elec- 
trodeless discharge. 


(3) Wheatstone’s bridge. 
(4) Set to purify commercial nitrogen. 


The measurement lamp was,built slightly different from Lord 
Rayleigh’s lamp (to facilitate exchanging the foil), see Fig. 1. The 


\ WOOO OC OC OOOO! 
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Fig. 1 


platinum foil had the following dimensions: width 1,4 mm, length 
about 26 cm and thickness 0,03 mm. The thin tube R may be placed 
on the foil by means of the ground joint 5, by appropriate tilting of 
the whole lamp. The joints S, and S, were smeared with Apieson grease. 
The volume of the lamp was determined to be 350 cm? and its dia- 
meter 5cm. The coil activating the gas in the discharge tube was 


as f 
` A `~ 
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iat the same time the self-inducting coil of the high frequency current 
ı circuit. As a condenser two Leyden bottles were used. Approximate 
‘estimates gave the value of the frequency of damped oscillations to 
be about 2 megacycles per second. 

Two resistances of 3M® each connected in parallel to the Ley- 
iden bottles served to remove the remnants of the charge and to sym- 
metrize the whole system. The maximum tension on the spark plug 
amounted to about 30 KV. A specially constructed interrupter (for 
a current of 50 A in the primary circuit of the inductor) powered 
by a synchronic motor gave two interruptions of the current per 
second. 

The nitrogen was taken from a steel flask. Water and carbon 
dioxide were removed in a tube filled with natrium lime; next a tube 
filled with pure copper heated to about 500°C removed the oxygen. 
P,O; absorbed remnants of water, and before entering the measuring 
lamp the nitrogen passed through a freezer cooled with liquid air. 


3. The results 


If we put R — electric resistance of the foil, Jį — current pas- 
sing through the foil while unscreened (in A), J, — dtto while seree- 
ned, p— pressure of gas in mm Hg, v— volume of the lamp = 
350 cm’, T — room temperature = 20°C, then the energy per particle 
of gas in the lamp (in eV) may be written as 


taR 1) 


Eey = 28p eV 


The following table gives the results of one series of measurements. 


No. current through inductor pmm J,A J, A RE eV/particle 


J 31 0,27 0,405 0,470 0,963 0,055 
2 31 0,188 0,46 0,504 0,995 0,06 
3 31 0,114 0,49 0,537 1,047 0,11 
4 31 6,8.1072 0,517 0,55 1,13 0,16 
5 31 4,07.10 2 0,528 0,562 1,235 VEG 
6 31 2,38.10—2 0,541 0,560 1,36 0,33 
7 31 1,2.1072 0,542 0,553 1,47 0,39 


Although the above results differ from those obtained by Lord 
Rayleigh, I consider that the differences are due to differences in the 
construction of the lamps and perhaps to differences in the discharge. 
On the other hand changes of energy depending on changes of pressure 
agree with Rayleigh’s results. 
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Very high values published by Rayleigh lead to the assumption 
that they cannot be attributed to nitrogen. Having obtained values 
of the same magnitude for oxygen and hydrogen, Rayleigh does not 
explain their source and writes. „It is difficult to see how the result 
can be explained by current conceptions”. 


4. Elecrodeless discharges 


Electric conditions in the lamp during the measurements were 
analyzed. These considerations are only approximate because collisions 
of electrons and ions with gas particles were not taken into conside- 
ration and because the structure of electric and magnetic fields was 
only approximately examined. The precision of results grows as pres- 
sure diminishes. The approximation may be sufficient to estimate the 
phenomenon because the great values of energy appeared when pres- 
sure was small and the mean free path was of the magnitude of the 
radius of the lamp or bigger. 

. The distribution of the electric and magnetic fields in an elec- 
trodeless discharge is simple in the case of a pure condenser or pure in- 
duction discharge1. The discharge in a eylindrical lamp on which 
a coil has been wound through which there passes a current of high 
frequency is of interest. It is a type of a ring electrodeless discharge. 
This discharge is at the same time of both the condenser and induction 
type. The simplest phenomenon occurs in the central plane of the 
coil vertical to its axis. The discharge is effected by the condenser 
field parallel to the axis and the induced field in a plane vertical to 
_ the axis. Both fields are in phase with each other. Owing to this, the 
direction of the resultant field -will not lie in a plane vertical to the 
axis aS would be the case if only the induced field existed. The di- 
rection of the resultant field depends on the ratio of the field strengths 
of both fields. The field strength of the condenser depends in a high 
degree on the extent of the surface of the conducting material of the 
coil. If the coil is very closely wound, it can be assumed that the elec- 
tric field is nearly homogeneous. We may say approximately that 
the condenser field F, of the coil will have the same value in its cen- 
tral part as if the whole potential difference U on the terminals of 


the coil were divided evenly over the whole length l of the coil, i. e., 
E,~U/l. 


1 The electrodeless condenser discharge is due to a varying electric field in 
a condenser; the induction discharge — to an electric vortex field which arises from 
a varying flux of magnetic induction. 
~ Cs 
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If we do not consider the Ohm resistance of the coil, which is 
comparatively small in comparison to its impedance (at frequencies 
of some megacycles per second), then we can assume that the electric 
current in the coil changes so as to make the EMF appearing in the 
coil equal and opposite to the potential difference applied. 

We may assume approximately that the potential difference per 
loop is U/n. This potential difference is evenly distributed over the 
whole length of the loop 2xR (R — radius). In optimum conditions 
the induced field strength H, (in external parts of the lamp) equals 
approximately U/2anR; in a distance r from the axis E,~U,/2anR?. 
Thus H,/E,~2ank?/lr; e. g. for 1=2R, and n=10. The field strength 
of the condenser field will be 10x times bigger than the induced 
field strength. This means that the condenser field strength will be 
higher by one order of magnitude than the induced field strength. 
It is obvious from the above that the influence of the condenser 
field will grow as the winding becomes closer and the radius larger 
(provided l> R). 


3. Possible explanation of Rayleigh’s results 


Rayleigh considered only the induced electric field in spite of 
the fact that the geometric conditions of his lamp facilitated the ap- 
pearance of the condenser field. Furthermore, he assumed that the 
nitrogen particles could not be activated in the nearest vinicity of 
the foil, which was placed in the axis of the coil, where the induced 
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=o 


Fig. 2 


field strength vanishes. But the nitrogen in the vinicity of the foil 
could be activated by rapid electrons accelerated in external parts 
of the tube, and deflected to the axis by the magnetic field of the 


- coil. The electrons themselves striking on the foil must have had a much 


greater influence on the amount of energy received by the foil. Ac- 
cording to calculations electrons — but not ions — could provide as 
much energy as that measured by Rayleigh. These energies are large 
for pressures under which ihe mean free path of the electron is equal 
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to, or greater than, the radius R of the tube. In Rayleigh’s lamp Æ 
was of the order of 1 cm, and the energy per particle rapidly rose 
from 0,1 mm Hg on. 

The introduction of the conducting foil did not affect the in- 
duced electric field, but it must have completely changed the con- 
denser field. If in the lamp without the foil the condenser field was 
almost parallel to the axis (except at both ends of the coil), after the 
introduction of the foil as an equipotential plane its direction was 
almost perpendicular to the axis. This indicates that in the sttong 
field between the foil and the coil an intensive discharge took place, 
supplying large amounts of energy to the foil. Approximate calcula- 
tions prove that in this effect under the conditions of Rayleigh’s 
work ions and electrons must be taken into account. This is in accord 
with R. E. Worley’s note (1948). It also follows that several other 
effects did not have sufficient mean strength. In the present work 
the energies obtained were smaller owing to different geometric and 
electric conditions in the lamp used (Fig. 1). 

The writer wishes to repeat his expression of gratitude to Pro- 
fessor J. Niklibore for suggesting the problem, for valuable discus- 
sions, and for his attentive care during the carrying out of this work. 
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A Note on the Franck-Conden Principle 


A. JABLONSKI 
Physics Department, Nicholas Copernicus University, Torun 
September, 9, 1951 


It appears that the quantum mechanical version of the Franck-Condon prin- 
ciple can be obtained by applying the „sudden approximation“ method to the eigen- 
functions describing the nuclear motion in a system consisting of electrons and 
nuclei, say, in a molecule. To show it one has to use the formalism given e. g. in 
Schiff’s Quantum Mechanics (New York 1949) in the chapter dealing with the me- 
thods for time-dependent problems. The conditions of applicability of this method 
to our particular problem are much the same as those for which the first Condon 
approximation is valid, the approximate separability of the electronic and nuclear 
motion being the most important. The eigenfunctions for nuclear motions have to 
be obtained from Schrödinger equations involving different potential functions 
corresponding to different electronic states of the system. 

Let an electronic transition occur in a system at t= 0. Since at t=0 the po- 
tential function for nuclear motion undergoes a sudden change, the , nuclear“ Ha- 
miltonian H, for t>0 will differ from the nuclear Hamiltonian H, for <0. In other 
words, the nuclear Hamiltonian H(t) changes discontinuously from H, to H, at t=0. 
Thus the „sudden approximation“ can be applied. From the equations 


Hyun =Entn t<0 3 (1) 
H, vm =Emvm t>0 (2) 
(the notations are those of Schiff) one obtains two complete orthonormal sets of 


nuclear eigenfunctions corresponding to two electronic states of the system. 
The general solutions of (1) and (2) are 


iEn a 
TERS] j r 
p=} anune h > t<0 ; (3) 
n 
iE, á 
sey 
Nomim A to. (4) 


Since the wave function must be a continuous function of the time also at 
t=0, (3) must be equal to (4) for t=0: 


2 antin= 3) bmm. (5) i 
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Thus 
bm asi) an fom undr, (6) 
n 


|bm |? being the probability discribution of the quantum states „m“ of nuclear mo- 
tion for t>0, provided the distibution for t<0 be given by aĝ (n,m = 0,1,...). 

If the system is in a particular state, say n, at t<0, (an=1, all others a’s equal 
to zero), (6) becomes 


bnm = [om Un dr. (7) 


The square of (7) 


(eam i= | fem un dz (8) 


is the probability of a change n-m of the state of the nuclear motion. In other 
words, (8) gives the relative intensity of the different n->m transitions associated 
with a definite electronic transition. Eq. (7) and (8) are identical with those resulting 
from the quantum mechanical version of the Franck-Condon principle}?. 


The Derivation of Vlasov’s Equation from Fock’s Equation 


JAN LOPUSZANSKI 


Institute of Theoretical Physics, Wrocław University, Wroclaw 


October 11, 1951 


The classical theory of plasma was given by A. A. Vlasov?; its basic equa- 
tion is 
of — — —; ra d ¢ # tA 
31 tE Val tPF Zp Vpl Te |K iae |) fap f(@, p')=0, A) 


where f=f(r,p,t) is the distribution function of electrons, p is the momentum, F the 
e2 
|a—a’ | T Ta, g) 
the mutual interaction of electrans is taken into account, but. changes of f due 
to ionization, recombination and scattering of electrons are disregarded. 
Equation (1) can be obtained from approximate equations of quantum me- 
chanics for identical particles. M. A. Mokulskij? obtained Vlasov’s equation by 
Hartree’s self-consistent field method by the limiting transition h>0. Mokulskij 
considers only a system of two electrons without taking into account either the 
exchange energy of two electrons with the same spin, or the influence of the ex- 
ternal field. It seems necessary to impose certain conditions on the function of 
potential energy to ensure correctness of the reasoning. These conditions are that 
the function be linear and that it does not depend on p. Mokulskij takes Terletskij’s 
function ê as the quantum correspondent of the classical distribution function: 


external force, K(|*—<’|)= Coulomb’s energy. In this equation 


2i ipx 
f = (2mh) * w*(a,t)g(p,t)e”. 


1 Condon E. U., Phys. Rev. 32, 858 (1928). 

2 Vlasov A. A., J. Exp. Theor. Phys., USSR, 8, 291 (1938). 
è Mokulskij M. A., ibid., 20, 688 (1950). 

4 Terletskij J. P., ibid., 7, 1290 (1937). 
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In this paper I endeavour to prove that Vlasov’s equation may also be obtained 
from the principles of quantum mechanics by another method. The point of issue 
is Fock’s equation’. The derivation is based on Dirac’s method of the density 
matrix? and on passing to classical mechanics by means of Bohr’s correspondence 
principle. A similar method was employed for instance in obtaining the Thomas- 
Fermi-Dirac equation %. In this method (i) exchange energy is considered, (ii) no 
conditions must be imposed on the potential energy, (iii) the influence of external 
fields is taken into account. 

Fock’s equation in the coordinate representation takes the following form 


n 


fc “| Hy | 2”) de” pria”) + f| P Bjj| x 
j=1 


PICCE ) dx” 97 ( a 


DEAN all’ 
5 t >. at ut 2 pj (a ) tt 
where (xBjr| z) = Zye ERA nae da’, 


77s ) dz” prls oi) 


on By signifies the energy of interaction, D Bjk the exchange energy, gz(a) the 
F 


i 
wave function of the electron, n the number of electrons in the system. The den- 
sity matrix in the same representation aN be written 


(æl e |2”) -Sae Plz) ¢ (3) 


assuming that each state pz is equally probable and that 
f(elelxaz =n 


takes place. By differentiating equation (3) partially with respect to time and sub- 


* 
stituting for - and ZPL the values from equation (2) we obtain 
h 3 A A P “st 
—=7 (2'|0|2") = (a'|Ke— eK |x”), (4) 
i ot 
k 
where K=H, a B— A is the energy operator; H, = — Oo y? + U(x), 
CHICA] we") yl h i ‘ p E a(x lolx”) 
B= - 3% = e SE ITT is the energy ofinteraction, and (x’|.A|x’’) ra, a) 


the exchange energy. If we treat K as the energy operator of the whois system, 
we may write the above equation simply as 


eee 


where Pee 2e Tar [k, elg (to. blo = = (ab — ba)) 


(5) 


1 Fock V., Z. Phys., 61, 126 (1930). 
2 Dirac P. A. M., Proc. Camb. Phil. Soc., 26, 376 (1930); Frenkel J. I., Wave 


Mechanics (in Russian), II, Leningrad 1934. 
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3 Gombás P., Die statistische Theorie des Atoms, Wien 1949. 
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To be able to pass to classical values we express the arbitrary dynamic va- 
riable Q in a mixed coordinate-momentum representation: 


ixp 


(lQlp) = f@1Qla’) aa’ e” 


and make the following assumptions: (i) we do not take into account the spin va- 
riable on the condition, however, that the exchange takes place only between elec- 
trons of parallel spins, (ii) the values p and x commute, (ili) n is very large, which 
involves (x|: |x) aa Rett ). Expounding the classical correspondent of the dy- 
namic variable Q into a Fourier’s integral and employing the correspondence prin- 


ciple, we obtain 


(6) 


ip(x/—x) 
I dx. i 
cene faoa * ™ ci 
Comparing equations (6) and (7), we see that 
ipx 
a (8) 
(2|Q|p) = Q(x, p) e 
(21Qlz) = (224? [Q(@,p)a (9) 
Equation (5) will transfer therefore by correspondence into the equation 
de (2, p 2 
LEP) + [E (a, p), ele, Pile = 0 (10) 


where [a,b], is the Poisson classical bracket expression, and ọ(%,p) represents a cer- 
tain classical state. Substituting the function @ for Q in equation (9) we obtain an 
interpretation of o(z,p) as the number of electrons in the volume (2zh)*. The dis- 
tribution function will therefore be expressed by 


f(x,p) = (2ah)~ e(a,p). (11) 
In order to obtain K(p,7)=H,(p,x) + B(p,x)—A(p,x) we have to calculate sepa- 
; 2 
rately each expression. Thus H(x,p) = Pa + U(x), 


2m 
1 
B =. —3 Ji 1 , KE} 
(x) = (2ahy Set f gr f dp’ e's’) (12a) 
Introducing spherical coordinates in A as integration variables we obtain 
2 / z 
DG ea E 
(x, p) Tabala dp’, (12b) 


where the fact that the exchange takes place only between electrons of parallel 
spin has been taken into account. 


From this as an equation for f we obtain 
of Dorie c Mle we 
Spal ele TENT I V N e fas aay f WP fas’) (13) 


2mh? J LART ETA Mi Lowa : 
+ 2% oy fp pp? p)dp A fa AA as 0, 


where F= — G u(x). 
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If in this equation we disregard the terms containg h?, we obtain Vlasov’s 
equation (1). If the integrals in these terms are large, or Vent is large for certain 
values of p. these expressions should not be disregarded. In this case Vlasov’s equa- 
tion takes the form of a differential equation. 

I render sincere thanks to Dr R. S. Ingarden for his help in this work. 


Electrostatic Neutron-Electron Interaction 


J. OLSZEWSKI 


Nicholas Copernicus University, Torun 


November 27, 1951 


The problem of the static interaction between neutron and electron is in- 
vestigated. This problem was investigated experimentally and theoretically, but 
the results of several authors! did not agree with each other. 

In particular Case and Borowitz & Kohn divided the „effective current“ 
of the neutron into two parts: 

ôju = bjut bi, (1) 
where the first part describes the current connected with the spin (and the ma- 
gnetic moment) of the neutron, while the second part is spin independent. For the 
evaluation of the neutron-electron interaction the above mentioned authors con- 
sidered only the spin independent part 6 i The first part Shy was neglected with 
the remark that the electric dipol moment associated with the magnetic moment 
of the neutron is negligeably small in the static limit (very slow neutrons). 

A direct computation showed, however, that by calculating the considered 
effect without explicit partition (1) into bin and by one obtains results incom- 
patible with the above mentioned authors. On the other hand, by taking account 
of the neglected current Sin: our results may be brought into agreement with the 
results of Case. 

The computation of the neutron-electron interaction is greatly simplified 
in the special case of central collisions. At the same time the assumption of central 
collisions seems very natural for the evaluation of the static potential. Namely, 
we can show quite generally that the frame of reference connected with the mass- 
centre of the neutron (before and after the central collision) is identical with the 
true electrostatic frame of reference (i. e., the frame in which the electromagnetic 


> 
field of the electron is the static Coulomb field: ra =0, Ay=ip=—ie/d4z|r|). In this 
frame of reference the only component of the neutron „current“ is ją- Both parts 
of 6j, may be represented as a power series of (AK,,)?, where AK, denotes the dif- 
ference between the initial and final values of the momentum-energy vector of the 


neutron. 
The zero-order term of oj! means an electric charge of the neutron, whic 
can be removed by renormalization. This term may be also removed by the regu- 
# ; 
1 Slotnick and Heitler, Phys. Rev., 75, 1645 (1949); Case, Phys. Rev., 
76, 1 (1949); Dancoff and Drell, Phys. Rev., 76, 205 (1949); Borowitz and Kohn, 


Phys. Rev., 76, 818 (1949). 
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larization of the mesonic and nucleonic masses, which can be showed by direct 
calculation. The next term, proportional to (AK ,,)*, describes a polarization charge 
induced by the Coulomb field of the electron. This term may be interpreted as the 
effect of a partial dissociation of the neutron into a virtual pair: proton — negative 
meson forming an electric dipol. As the Fourier component of the Coulomb field 
is proportional to (AK, pas the transition to the limit AK -> 0? yields a finite result 
for the effective potential of interaction. 

Similarly, the zero-order term in the power expansion of ôj vanishes since 
there is no electric moment associated with the spin of a neutron at rest. On the 
other hand, the terms of the order (AK „)? do not vanish (owing to the Fourier com- 
ponent (AK ee of the Coulomb field) and yield a contribution in the limit 4K Pie a 


V 
AEV] 


Fig. 1 


Thus, besides ee spin ECT PER vonsidered i Case and Borowitz 
& Kohn we get also spin dependent terms from sj! 

The total effect may be represented with the aid of the depth of a square 
radi i width ro pori gives the same Papes aen (ro its classical electron 
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The technique used in my paper is that of Dyson with the use of regulators 
A detailed account will be published in this Journal. 

I should like to express ,my deepest thanks to Professor J. Rayski for the 
suggestion of this investigation as well as for much valuable criticism. I am also 
much indepted to the Komisja Popierania Twórczości Naukowej i Artystycznej 
przy Prezydium Rady Ministrów for a grant-in-aid. 

Remark. In the mean-time appeared a note of Foldy! in which the spin 
dependent part of the neutron-electron interaction is derived from the observed 
magnetic moment of the neutron. Thus, it is the same correction as our ój! obtai- 
ned in a different way. Foldy does not present details of his calculations. 


A Multiple Shielding Counter 
J. JOZWIAK and W. MOŚCICKI 


Institute of Experimental Physics of the University of Poznań 


March 13, 1952 


Seeger as well as Jürgensen and Weiss ? describe a Multiple Shielding Counter 
consisting of two coaxial cylinders and 6 anode wires enclosed in the space bet- 
ween the cylinders. Independently of the above mentioned research a counter al- 
most identical in principle was built and partly investigated in our Institute? as 
an anticoincidence shielding of internal carbon-dioxide filled counter for research 
carried by one of us (W. M.) as well as an anticoincidence shielding for neutron 
counters. 

Our counter consists of two coaxial cylinders (diameter d=45,5 mm, and 
D=75 mm), it contains n,=8 anode wires (2r=0,1 mm) and n,=8 cathode wires 
(2r=0,25 mm), -placed on the circumference of a circle (2r=65 mm) around the 
chief independent counter coaxial with the cylinders. The length of the counter 
is L=40 cm. Fig. 1 shows schematically the distribution of the electrostatic po- 
tential for each anode wire. 

The investigation of each counter wire separately revealed that the diference 
of threshold voltages of all wires lie wihin the region of 20 V, the plateaus extend 
over 200 V. ; 

The springing of the discharge between neighbouring counter-wires when 
operated over a common resistance, the phenomenon described by Jürgensen and 
Weiss, has never been observed. The number of impulses recorded in two neigh- 
bouring counters equals — within the limits of statistical error — the sum of the 
numbers recorded for'each counter separately. Very rare large sharp impulses of 
various amplitudes could be observed by means of a cathode rays oscilloscope. 


1 Foldy L. L., Phys. Rev., 76, 1 (1951). 

2 Seeger A., Z. Naturforsch., 6a, 688 (1951); Jirgensen H. and Weiss 
H. M., ibid., p. 694. = 

3 The idea of the Multiple Shielding Counter was conceived in autumm of 
1950 (J. J.). Difficulties in providing indispensable materials delayed its final con- 
struction till autumn 1951. However, our preliminary results proved soon promising 
enough to allow to put the problem of a systematic investigation of the properties 
of. M. S. C. in the official plan of research to be made in 1951/52 in our Institute. 
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We are inclined to ascribe them to coincidences of discharge in two or more coun- 
ters. Jirgensen and Weiss show that the efficiency of M. S. C. bombarded with 
fast g particles varies from 96°/, (for electrons moving in the plane A) to 99,8°/° 
with an average efficiency of 98°/,. They did not find the number of impulses to 
depend on the azimuth of incidence when M. S. C. were illuminated with a colli- 
mated beam of y rays. These properties do not seem to be of sm importence nor 


a ee ee 


decisive for the use of this maces as an anticoincidence shield considering the disg 
harge mechanism initiated by y rays and electrons. : 
We attempt to investigate the efficiency of the shielding function of our 
M. S. C. against the hard ionising component of cosmic rays. 
rf The comparison of our results with those of Jirgensen and Weiss shows thad 
the cathode wires yield certain advantages, though their introduction is not free 
from some technical difficulties. The cathode wires diminish the region of low sen- 
sibility and yield a clearer delimitation of the regions of the neighbouring counters. 
The lack of efficiency of M. S. C. in plane A (Fig. 1) can be eliminated ty 
ducing an odd number of anode wires. Details on the construction of a counter 
this type with an odd number (n=7 or 9) of anodes and es as men as 
on the prerana of its properties, sol Be published. 
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ERRATA 

A Note on the Franck-Condon Principle, A. Jabłoński, Acta Phys. Polon., 
11, 195 (1951). | 

_ Equation (4), p. 195, should read 
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